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ABSTRACT 


Anomalous  atmospheric  conditions  at  Western  Test  Range  were  measured 
bv  aircraft.  Tha  invarainn  layer  was  found  to  possess  upward  wsui*p  as  it 
approached  the  coast,  as  well  as  dilation,  lensing,  sandwiching,  and  waving. 
Flow  geometry,  when  associated  with  topographic  obstacles,  was  found  to 
be  markedly  similar  to  that  produced  on  laboratory  models.  Froude  number 
concepts  classified  the  anomalies  accurately.  Refraction  corrections  pro¬ 
duced  by  three  possible  data  sources  were  compared  to  assess  the  effects 
of  the  above  described  anonaalies.  When  data  from  full  cross  sections  along 
the  sighting  path  were  employed,  along  with  steering  equations,  elevation 
and  bending  corrections  were  considerably  greater  (2x)  than  those  produced 
by  single  point  radiosondes  or  aircraft  spiral  measurements.  Moreover, 
corrections  produced  by  spiral  data  increased  to  seaward  where  the  inver¬ 
sion  became  sharply  defined.  The  validity  of  full  cross  sections  defining 
the  propagation  path  was  best  established  by  examining  ray  paths  for  antenna 
elevation  angles  existing  at  the  time  of  signal  loss  when  a  target  simulation 
aircraft  encountered  a  radio  hole.  Ideally,  a  trapped  ray  could  be  expected 
to  bend  back  to  earth  at  the  aircrafts'  position  in  time  amd  space  when  any 
one  of  three  tracking  radars  lost  lock  on  the  target.  A  cross-section  corr¬ 
ected  propagation  path  defined  this  target  position  better  thaui  any  other 
means  of  defining  the  path.  Strength  and  ^stinctive  geometry  of  the  inver¬ 
sion  at  the  position  of  loss  are  displayed.  Mathematics  employed  for  real¬ 
time  corrections  at  any  angle  (including  negative)  are  presented,  as  are 
mathematics  for  steering  corrections. 
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INTROnilCTION 


The  location  of  this  investigation  is  Western  Test  Range.  The  purposes  of 
the  investigation  were: 

(1)  Define  the  hazards  (anomalies)  of  the  atmosphere  with  respect 
to  their  bearing  on  tracking  problems. 

(2)  Define  tracking  accuracies  at  lower  limits  (elevation  angles  zero 
to  5  degrees). 

(?)  Assess  currently  employed  refraction  correction  techniques 
answer  the  question  "How  are  we  doing?" 

Structure  of  the  Report 

Citing  the  structure  of  this  report  will  more  easily  allow  one  to  identify  the 
logic  and  interdependence  of  topics  to  be  presented.  Appendix  B  contains  a 
review  of  past  work  on  the  relation  between  the  north  Pacific  inversion  and 
radio  propagation.  The  mechanisms  of  the  inversion  are  explained,  the  fac¬ 
tors  controlling  its  activity  (geometric  changes  and  strength  changes),  are 
discussed,  and  field  measurements  of  its  effect  on  radar  are  reviewed. 

Section  II,  Description  of  Atmospheric  Anomalie.s,  explains  the  dominant 
types  of  anomalies  found  at  Western  Test  Range.  Their  causes  and  spatial 
distribution  are  explained.  Section  III,  Comparison  of  Refraction  Correction 
Input  Data,  compares  amounts  of  correction  applied  to  a  given  tracking  event 
as  produced  by  three  possible  data  sources,  balloon  soundings  (Raob);  air¬ 
craft  spirals  and  fully  recorded  refraction  cross  sections  along  the  sighting 
trajectory.  Since  these  correction  outputs  are  relative  to  each  other.  Section 
IV  is  devoted  to  providing  an  absolute  reference  for  deciding  which  of  the  three 
is  closest  to  reality  and  how  close. 

Since  we  were  unable  to  get  the  Western  Test  Range  supplied  REFRAC  com¬ 
puter  program  to  run  through  the  computer,  we  re-prograunmed  the  correc¬ 
tion  calculations  using  the  original  mathematics  to  maintain  consistency.  A 
review  of  assumptions  and  possible  errors  is  presented.  Since  this  program 
is  unable  to  handle  the  anistropic  atmosphere  existing  at  Western  Test  Range, 
a  ray  steering  program  was  developed.  In  connection  with  this  effort,  a  real¬ 
time  correction  technique,  good  at  any  angle  including  negative  angles,  was 
produced.  Section  V  presents  the  structure  of  these  programs  and  culminates 
in  the  handling  of  steering  problems.  Appendix  A  contains  supplementary 
.mathematics  for  Section  V.  Section  VI  contains  conclusions  and  recommenda¬ 
tions. 
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SECTION  II 


DESCRIPTION  OF  ATMOSPHERIC  ANOMALIES 


In  defining  the  atmospheric  anomalies  existing  at  Western  Test  Range 
(including  the  San  Nicolas  Island  and  Point  Mugu  tracking  stations),  we 
wish  first  to  introduce  work  of  Robert  Long  (1953,  1954,  1955)  in  flow 
dynamics,  which  discusses  understanding  the  flow  features  pertinent  to 
the  anomalies  fovind  at  Western  Test  Range.  In  the  present  instamce  we 
use  the  term  "anomalies"  to  define  the  perturbations  in  the  atmosphere 
wdiich  axe  characteristic  of  its  departure  from  either  a  well-mixed  stand¬ 
ard  atmosphere,  or  from  the  often  genersdized  stratified  atmosphere.  A 
stratified  atmosphere  can  be  described  as  an  inversion  layer  at  about  the 
2,  000  ft  level  with  moist  marine  air  below,  and  subsiding  dry  air  aloft. 

As  a  generalized  concept,  this  is  correct.  As  a  definition  of  the  specific 
medium  through  which  missiles  are  tracked  at  Western  Test  Range,  it  is 
not  correct  when  one  considers  low-angle  tracking  with  very  accurate 
radar.  Early  low-angle  acquisition  of  incoming  missiles  could  be  a  pract¬ 
ical  reason  for  wanting  to  know  more  about  anomalous  propagation,  as 
could  the  tracking  of  low-altitude  aircraft.  Before  further  description  of 
atmospheric  anomalies  is  discussed,  the  concept  of  Froude  numbers  and 
flow  characteristics  should  be  grasped. 

The  Hydraulic  Jump 

When  an  obstacle  is  placed  in  the  path  of  moving  fluid,  a  wave  is  formed. 

The  geometry  of  this  wave  is  dependent  on  the  density  of  the  fluid,  the 
height  and  geometry  of  the  obstacle,  and  the  velocity  of  the  fluid.  These 
relationships  are  explained  as  follows  (.Skalak,  1967). 

In  open  channel  flow,  the  counterpart  of  the  Mach  number  is  the  Froude 
number  F,  defined  as  the  ratio  of  the  flow'  velocity  V  to  an  elementary  grav¬ 
ity  wave  velocity,  c  =  / gd,  where  d  is  the  depth  of  flow.  For  Froude  num¬ 
ber  F  =  V// gd  less  than  unity  the  flow  is  called  subcritical  and  is  ir,  many 
ways  analogous  to  subsonic  flow. 

Supercritical  is  defined  as  open  channel  with  Froude  number  F  >  1.  It  is 
analogous  to  supersonic  flow.  Flow  in  a  tranquil  river  is  usually  subcritical, 
but  rapids  are  super criticaL  Lopez  and  Howell  (1967)  wrote  of  Froude  num¬ 
ber  application  to  katabatic  flow  in  the  Andes.  The  same  is  applicable  to 
flow  over  Tranquillon  Peak,  the  Vandenberg  tracking  site;  and  San  Nicolas 
Isltuids,  sdl  of  which  are  units  of  the  Western  Test  Range  trilateration  sys¬ 
tem. 
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Example  of  Air  Plow  Across  a  Ridge  with 


Figure  1 
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to  designate  initial  flow  velocity,  gravity,  amd  wave  height,  Lopez  and 
Kowell  (ibid)  explain  airflow:  When  cool  air  is  considered  instead  of  water, 
gravity  acts  on  the  density  differences  between  the  air  strata  in  accelera¬ 
tion  the  flow,  and  g  in  Eg  (1)  is  replaced,  at  least  for  relatively  small 
vertical  displacements,  by 

y  =  (M/B)  .  g  (2) 

where  B  is  no  longer  constant  but  dependent  on  the  relative  potential  temp¬ 
erature  difference  A0.  The  critical  spill  velocity  Eq  (1)  then  becomes 

V-  =  (yh  (3) 

o  o 

Let  us  compare  two  cross  sections  0  and  1  (Fig.  1),  assuming  a  frictionless 
streamline  flow.  If  H  is  the  height  of  the  pass  over  the  valley  floor,  apply¬ 
ing  Bernoulli's  equation  to  cross  sections  0  and  1: 

y(H  4h^  -  h,)  =  1/2  (v^  -  v^^)  (4) 

The  equation  of  continuity  requires  that 

Discharge  =  Q  =  Vjhj  “  (5) 


I 
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Combining  Eq  (4)  and  (5)  with  the  critical  velocity  of  Eq  (3)  we  arrive  at 

(h  /h,)^  -  (2H  /h  +  3)(h  /hj  +  2  =  0  (6) 

o  1  o  o  O  1 

The  Froude  number,  F.  -  v^/v  expresses  the  degree  to  which  the  flow  is 
supercritlcaL  The  ratio  (h^/h^  is  always  less  than  1  and  equals  the  shrink 
ratio  between  the  depth  of  the  supercritical  shooting  flow  at  the  base  of  the 
pass  and  depth  of  the  overflowing  air. 

As  the  valley  floor  of  Fig.  1  flattens,  the  supercritical  velocity  of  the  flow 
is  rapidly  diminished  by  friction  and  a  "tranquil"  subcritical  flow  region  is 
approached.  In  order  to  satisfy  continuity  requiren^ents,  the  depth  of  the 
flow  must  increase  from  to  a  certain  h^.  If  this  happens,  waves  can 
travel  upstream  in  the  subcritical  area  near  h^  but  cannot  do  so  in  the  super¬ 
critical  area  above  h^.  The  wave  front  will  therefore  steepen  until  it  breaks 
down  somewhere  between  cross  sections  1  and  2.  In  this  way,  a  hydraulic 
jump  is  formed  which  tends  to  remain  stationary.  A  simple  illustration  of 
the  hydraulic  jump  may  be  seen  by  rxmning  tapwater  into  an  empty  basin  (Fig.  2). 


4 


"'lllmgl'illi'.-  lb:*  .  .1.1. in...  'Jill  . . .  . . I  •  ..  .  ■. .. . . 


f 


I 


p 

[£- 


& 


I 


V.Ticrc  the  streani  hits  the  basin;  a  thin  sheet  of  fast-raovlng  water  with 
supercritical  velocity  is  surroxmded  by  a  ring  of  deeper  turbulent  water. 

A  pencil  point  set  in  the  thin  sheet  will  creat  a  "shock  wave"  of  the  familiar 
V-shape  which  clearly  does  not  propagate  upstream.  This  jump  is  often 
referred  to  as  the  "pressure  jump"  when  it  occurs  in  the  air. 

Illustration  of  the  above  relationship  of  these  parameters  was  supplied  by 
Long  <1953,  1954,  1955)  and  in  motion  picture  experiments  where  he  showed 
the  flow  geometry  of  waves  in  a  two -liquid  medium  of  two  densities  as  the 
fluids  flowed  over  the  symmetrical  obstacle.  Long  constructed  an  elongated 
water  tank  similar  to  tropical  fish  tanks  kept  by  hobbyists,  in  which  he  irtioto- 
graphed  a  moving  obstacle  and  the  wave  caused  by  the  obstacle.  In  this  case, 
the  obstacle  moved  instead  of  the  fluid,  the  results  being  the  same.  The 
camera  moved  with  the  obstacle,  photographing  the  wave  as  it  proceeded. 

The  geometry  of  the  obstacle  was  varied,  that  is,  the  height-length  relation¬ 
ship,  as  was  the  velocity  of  the  obstacle  through  the  fluid.  Figures  3,  4,  and 
5  extracted  from  16  mm  movies  of  Long's  work,  show  how  the  shape  of  tije 
wave  is  changed  with  increasing  Froude  numbers  in  a  two-density  stratified 
system.  Figure  3  shows  a  wave  train  in  the  lee  of  the  obstacle  which  is  gen¬ 
erated  at  low  Froude  numbers  in  the  denser  layer  (low  fltiid  velocities).  Fig. 

4  shows  the  convergence  of  the  lower  medium  on  the  leeward  side  of  the  ob¬ 
stacle  as  the  velocity  increases;  and  shows  a  hydraulic  jump  has  occurred  on 
the  leeward  side  of  the  obstacle.  In  this  latter  instance,  the  flow  of  the  fluid 
is  faster  than  the  travel  time  of  any  wave  that  might  be  propagated  on  the 
leeward  side  of  the  obstacle  and  allowed  to  travel  upstream.  Since  this  energy 
cannot  pass  upstream,  it  collects  at  a  point  downstream  where  the  velocity 
drops  sufficiently  to  allow  wave  generation.  At  this  point,  we  have  a  hydraulic 
jump.  With  still  higher  velocities  a  smooth  symmetrical  wave  forms  in  the 
stream  over  the  obstacle.  Figure  5.  This  can  be  seen  in  nature,  for  example, 
as  a  standing  wave  manifested  by  a  symmetrical  lenticular  cloud  over  the 
Sierra  Nevada.  Referring  still  to  observance  in  nature,  we  will  see  in  the 
following  work  that  the  nature  of  the  marine  layer  and  the  less  dense  air  aloft 
at  the  Pacific  Missile  Range  behave  much  like  the  idealized  laboratory  models 
seen  in  the  photographs.  The  marine  layer  represents  the  denser  medium 
shown  in  the  photographs.  The  place  of  the  obstacle  is  taken  by  the  topograph¬ 
ical  land  mass  of  Tranquillon  Peak  or  San  Nicolas  Island.  The  inversion 
interface  is  the  line  of  demarcation  between  the  two  fluids  in  the  idealized 
laboratory  presentation.  We  will  see  from  direct  aircraft  measurements  that 
the  nature  of  the  flow  in  the  vicinity  of  the  tracking  sites  at  Western  Test  Range 
is  identical  to  that  of  the  laboratory  prediction. 

Aircraft  Measurements 


Aircraft  missions  were  designed  to  obtain  cross-sections  along  the  tracking  i 

line  of  site  between  Tranquillon  Peak  and  the  Minuteman  Missile  pads  20  1 

miles  to  the  north  of  Tranquillon  Peak;  also,  a  cross-section  was  taken  in  the  \ 
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FIGURE  3.  Longs  loborotory  photo  of  low  Froudo  nufnbw. 


FIGURE  4.  Longs  loboratory  photo  of  modsrtits  Frouds  nuntosr. 


FIGURE  5.  Long's  laboratory  photo  of  high  Frouds  number. 
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direction  of  westerly  missile  trajectories;  and  a  cross-section  wais  taken 
along  the  tracking  line  of  sight  for  north- south  orbiting  vehicle":  Similarly, 
flights  in  the  vicinity  of  Pt.  Mugu>  and  San  Nicolas  Island  were  flown.  In 
order  to  obtain  these  cross -sectionn.  the  C-lSl  airplane  equipped  with  ana¬ 


logue/digital  recording  refractometer  (cycling  time  =  0.  1  sec)  radio-altimeter, 
and  other  meteorological  sensing  devices  was  flown  at  successively  increasing 
500-ft  levels  »long  the  line -of- sight.  Variations  in  the  atmosphere  observed 
as  waves,  len&<.8,  sandwiclUng  of  layers,  etc. .  show'ed  up  in  the  variations 
of  the  refractometer  records.  The  inversion  interface  was  located,  and  the 
aircraft  maneuvered  through  the  inversion  repeatedly  (successive  dives  end 
climbs)  such  that  an  inversion  penetration  occurred  every  1. 200  ft  or  so. 

The  waves  that  were  measured  were  on  the  order  of  6, 000  ft  in  wave-length 
and  often  much  greater.  Therefore,  about  six  penetrations  per  wave  were 
obtained,  each  penetration  marking  the  altitude  of  the  inversion  interface  at 
that  particular  point.  Moreover,  the  aircraft  was  flown  horizontally  through 
a  wave  tradn  and  the  wave  fluctuations  appeared  on  the  refractometer  trace 
as  continuously  changing  refraction  values  for  a  constant  altitude.  With  these 
inputs,  it  was  possible  to  reconstruct  the  structure  of  the  atmosphere  along 
a  given  cross-section.  These  reconstructions  were  made  in  the  laboratory, 
and  the  theoretical  paths  of  radar  rays  (bundles,  waves)  were  traced  through 
these  disturbed  cross-sections  by  computer.  The  laboratory  ray  tracings  are 
intended  to  reveal  the  hazards  of  the  atmosphere  in  quantitative  form  (Sections 
II  and  III),  and  indeed  must  be  consistent  with  actual  signal  losses  experienced 
in  tracking  the  C-band  beacon  of  the  C-131  aircraft. 


Case  Studies 


The  following  case  studies,  in  the  vicinity  of  Tranquillon  Peak  and  San  Nicolas 
Island,  and  Pt.  Mugu,  show  the  wave  form  imposed  on  the  atmosph'-re  by 
thermal  and  topographic  disturbances  as  defined  by  our  airborne  measurements. 


Case  1-4  October  1968 


Figure  6  represents  a  cross-section  between  Tranquillon  Peak  and  the  nor¬ 
thern  Minuteman  Missile  pads.  This  section  is  a  good  example  of  a  low- 
velocity  wave  train  (low  Proude  number.  Figure  3)  triggered  by  hills  just  to 
the  northwest  of  the  missile  pads,  that  is,  the  wave  train  originates  to  the 
left  (about  3-4  miles)  of  the  cross-section  shown.  Surface  wind  velocity  for 
October  4  was  8  knots.  The  shaded  layer  represents  the  inversion  interface, 
that  is  the  strong  gradient  zone  between  the  marine  layer  and  the  dry  air  above. 
This  is  where  the  strong  refractive  bending  takes  place  and  is  critical  to  our 
problem.  The  significance  of  this  section  is  that  it  not  only  defines  the  hazards 
of  the  Tranquillon  Peak  Minuteman  line-of-sight  but  also  closely  approximates 
the  look-angle  for  System  IV  (within  +  15®). 
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Case  2  -  Modified  Section 

More  pertinent  to  the  northerly  look-angle  irom  f  ranquiiion  Peak  ia  F igiu-e 
7.  which  was  modified  from  a  previous  diffusion  study  (Smith,  1964).  The 
modification  cunuiaieu  only  in  changing  the  measured  parameter  potential 
temperature  to  refractive  index.  Since  absolute  value  of  refraction  needs 
only  to  be  approximated,  the  relative  distribution  of  refractivity  is  the  im¬ 
portant  property.  This  section  was  used  to  emphasize,  through  use  of  in¬ 
dependent  sources,  that  the  perturbations  desertbed  are  read  to  all  field 
workers  and  can  be  found  most  any  time  one  wishes  to  make  the  effort  of 
measuring  them:  see  also  Cermak  (1966).  who  shows  the  same  effects  in  a 
1:12,  000  topographic  model.  This  (Smith)  section  is  particularly  tiseful. 

The  strong  inversion  on  this  day  happened  to  be  at  an  altitude  where  the 
marine  to  dry  air  interface  engulfed  the  Tranquillon  Peak  radar.  This  ex¬ 
ample  shows  a  cross-section  with  stronger  wind  speeds  than  Case  1,  and 
is  a  rather  typical  case  where  the  inversion  layer  is  warped  up  the  mountain 
on  one  side  and  strongly  converges  in  a  steep  descent  on  the  lee  side  as  in 
Long's  laboratory  illustration  (see  Figure  4,  moderate  Froude  number). 

Not  only  ia  there  attenuation  of  ttie  gradient  layer  on  the  downwind  side  of 
Tranquillon  Peak,  but  the  hydraulic  jump  and  resultant  mixing  are  also  evi¬ 
dent  in  the  spreading  of  the  refractopleths.  Figure  8  shows  a  photograph  of 
this  precise  Froude  condition.  The  flow  is  brou^t  out  sharply  by  the  cloud 
formation  which  is  following  precisely  the  refractopleths  shown  in  Figure  7. 

In  addition,  Figure  9  shows  where  the  mixing  due  to  the  hydraulic  jump  man¬ 
ages  to  destroy  the  inversion  layer  locally.  The  photograph  shows  a  high- 
altitude  view  of  Point  Arguello  and  Sudden  Ranch  coastal  area.  The  clear 
area  to  the  right  in  the  photo  is  the  area  of  mechanical  mixing  where  the 
winds  are  rounding  Point  Arguello.  The  dry  air  above  the  inversion  is  mixed 
into  the  marine  layer,  dissipating  the  cloud  droplets.  The  sharp  inversion 
interface  normally  found  between  the  marine  layer  and  the  dry  air  above  is 
destroyed  here.  The  cloud  erosion  in  the  lower  left  of  the  picture  shows 
destruction  due  to  hydraulic  jump. 

Figure  10  is  an  aircraft  measurement  of  the  inversion  layer  along  the  south¬ 
erly  look  angle.  Note  the  washing  out  of  the  inversion  in  the  right  side  of 
the  diagram  due  to  the  mixing  in  the  lee  of  Tranquillon  Peak  as  previously 
explained.  As  one  proceeds  to  seaward  some  20  kilometers  the  inversion 
is  reinstated  in  all  its  strength.  At  40  kilometers  the  inversion  is  sufficiently 
strong  to  produce  trapping.  We  will  consider  trapping  later  in  Section  IV, 
on  radio  hole  missions. 

Case  3-30  March  1967 

Figure  11  shows  the  westerly  trajectory  from  Tranquillon  Peak  to  seaward 
with  the  perturbation  caused  by  Tranquillon  Peak  showing  very  clearly.  Note 
in  the  diagram  that  the  dashed  set  of  lines  shows  a  section  taken  45  minutes 
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FIGURE  7.  North  section  from  Tronquillon  Peak 
showing  moderate  Froude  nuntber. 


FIGURE  8.  Dynamic  flow  oyer  Tranquillon  Peak 
revealed  by  clouds,  corresponding  to 
moderate  Froude  number.  (See  FIGURE  7) 
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FIGURE  9.  High-altitude  aircraft  photo  showing 
mixing  zone  to  west  and  south  of 

Tronquillon  Peak— Point  Argue! lo. 
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tude  of  the  fluctuation  to  remain  constant  >^le  the  wave  nodes  shift  in  spatial 
position.  The  sijpiificance  of  this  is  that  from  the  standpoint  of  refractive 
bending  the  magnitude  of  error  in  either  case  is  likely  to  be  much  the  same. 

As  we  follow  the  inversion  to  seaward,  it  stabilizes  some  20  or  30  miles  out. 

By  "stabilizes"  we  mean  that  it  becomes  strongly  stratified  at  one  altitude  with 
a  very  pronounced  gradient  layer  (as  opposed  to  multi  layering,  leasing,  or 
complete  mixing).  These  distinctions  are  not  very  discernible  to  single  point 
balloon  soundings  on  the  beach  (Vandenberg  or  Mugu). 

Case  4-24  July  1967 

Case  4  is  a  measurement  of  the  atmosphere  over  San  Nicolas  Island  on  24  July 
1967.  Figure  12  is  an  east-west  section  along  the  axis  of  the  island  and  along 
the  flow  of  the  wind  measured  by  aircraft.  Notice  how  a  "bow  wave"  forms  at 
the  windward  end  of  the  island.  The  turbulence  effects  induced  by  the  topo¬ 
graphy  are  evident  as  a  mild  hydraulic  jump  in  the  right  half  of  the  diagram. 
Strong  convergence  of  the  gradient  layer  over  the  land  mass  again  is  identifi¬ 
able  for  the  inversion  interface  at  mid  diagram.  This  repeats  the  moderate 
Froude  number  phenomena  seen  in  Figure  7  for  Traiiquillon  Peak.  Con .  er- 
gence  of  the  inversion  gradient  in  this  manner  is  extremely  significant  when 
examining  refractive  bending,  as  is  the  general  tilting  and  warping  and  waving 
of  the  inversion  interface.  Generally,  the  inversion  was  found  to  be  300  to 
400  ft  lower  on  the  north  side  of  the  island  than  on  the  south  side  of  the  island. 

Additional  clarity  of  the  flow  pattern  around  the  island  is  seen  in  Figure  13,  as 
viewed  from  a  high-flying  aircraft.  A  wave  train  progresses  out  from  the 
"nose"  of  the  island.  It  is  especially  evident  on  the  near  camera  side.  Figure 
14  shows  sections  across  these  waves.  There  is  a  doming  of  the  air  stream 
over  the  o-nter  of  the  island,  with  two  sister  perturbations  (the  V-wake)  to 
either  side. 

Summary 

To  summarize  the  anomaly  high  points  found  in  these  case  studies,  the  follow¬ 
ing  generalizations  are  made. 

1.  With  prevailing  northwest  wind,  the  inversion  immediately  to  the 
south  of  Tranquillon  Peak  is  likely  to  be  destroyed  by  the  hydraulic -jump  effect 
associated  with  sufficiently  high  Froude  numbers.  This  will  often  be  noticed 
as  an  erosion  of  the  cloud  deck  at  inversion  height.  This  circumstance  should 
allow  excellent  target  tracking.  The  well-mixed  atmosphere  here  will  approx¬ 
imate  a  standard  atmosphere.  The  mixed  condition  should  exist  some  20  miles 
to  seaward. 


2.  Generally,  the  area  immediately  to  the  west  of  Tranquillon  Peak 
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24  JUIX  196. 
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FIGURE  13.  High-altitude  view  of 


island  with  bow  wake-like 
structure  seen  in  the  stratus 

deck. 
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DISTANCE.  NAUTICAL  MILES 


will  have  low-level  mixing  due  to  eddying  around  the  land  mass  itself,  in  which 
case  one  will  find  a  clear  sky  area  along  the  beach  where  the  has  eroded 

the  clovid  deck  {present  at  inversion  height).  At  higher  altitude  one  will  find  a 
set  of  standing  waves  from  the  disturbance  of  tite  peak  itself.  Occasionally, 
serious  tracking  errors  caxi  be  produced  by  strengthening  of  the  gradient  layer 
a  short  distance  to  seaward  ('^  10  -  20  nm). 

3.  In  the  northerly  sections  wave  tr  ains  across  the  Santa  Inez  Valley 
can  be  expected.  The  Inversion  may  also  have  a  wturped  geometry,  which 
.seemingly  climbs  upward  to  meet  Tranquillon  Peak.  This  warping  can  actually 
be  seen  by  a  ground  observer  on  the  north  hills  when  the  low  morning  sun 
scatters  off  the  haze  layer.  The  trapping  and  bending  significance  of  a  warped 
or  tilted  inversion  layer  will  be  discussed  in  the  section  on  radio  hole  missions. 

4.  A  slight  elevation  difference  in  the  inversion  will  exist  between  the 
north  and  south  sides  of  the  San  Nicholas  Island.  The  gradient  layer  will  bulge 
over  the  island  on  the  order  of  1,  000  ft.  A  wave  structure  will  propagate  out 
to  the  north  and  south  sides  of  the  island.  There  is  strong  convergence  of  the 
gradient  layer  eastward  of  the  mid  point  of  the  island,  and  a  hydraulic  jump 
induced  mixing  Of  the  inversion  layer  downwind  of  this. 

5.  The  inversion  at  Point  Mugu  may  be  washed  out  from  eddy  mixing 
associated  with  shear  where  the  sea  breeze  front  meets  the  north  westerlies. 

It  can  also  be  multi  layered  or  lensing  tmder  moderately  weak  stability  con¬ 
ditions  (to  be  seen  in  Figure  16,  Section  IV,  8  June  example).  The  seaward 
strength  of  the  inversion  may  be  experienced  close  along  shore  at  Point  Mugu 
when  strong  stability  conditions  are  present. 
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COMPARISON  OF  REFRACTION  CORRECTION  INPUT  DATA 


Having  established  the  complex  nature  of  the  atmosphere  at  Western  Test 
Range  (Section  I),  quantitative  estimates  should  be  made  of  the  differences 
in  various  tracking  errors  as  they  depend  on  the  source  of  their  calculations 
inputs.  Comparisons  made  by  Gardner  (unpublished)  show  there  is  little 
difference  in  corrections  to  be  made  as  a  function  of  formulae  used.  The 
reason,  of  course,  is  that  all  calculations  are  based  on  Snell's  Law.  It  will 
be  apparent  in  the  section,  however,  that  one  cannot  afford  to  ignore  the 
validity  of  the  inputs.  One  faces  the  question  "How  realistic  are  the  measure¬ 
ments  being  fed  into  the  refraction  correction  calculations?" 

In  order  to  keep  the  comparisons  valid  we  must  use  common  radar  elevation 
angles  (1°),  common  target  altitudes,  and  approximately  the  same  points  in 
time  for  a  given  day.  Tables  1-6  represent  measurements  within  an  hour  of 
each  other.  Altitudes  are  limited  by  the  spiral  height  of  the  aircraft.  Cal¬ 
culations  are  compared  at  this  minimum  height.  The  important  points  of 
conformance  are  that  targets  be  above  the  anomalous  section  of  the  atmos¬ 
phere.  In  all  cases  it  is  the  relative  difference  in  tracking  errors  for  a  given 
day  that  we  wish  to  isolate.  The  general  magnitude  of  error  is  of  secondary 
interest. 

The  Table  1  format  displays  range  error,  elevation  angle  error,  and  bending 
error,  as  calculated  for  a  12,  000  ft  elevation  refraction  profile  measured  on 
4  October  1968.  The  primary  comparisoii  is  for  the  line  of  sight  between 
Tranquillon  Pe2dt  and  the  Minuteman  pads  located  about  20  miles  north  of  the 
Peak.  The  cross-section  measurements  produce  bending  corrections  about 
30  percent  greater  than  the  Vandenberg  Raobs  produce  (7.  2  versus  5.  6  mils). 

The  65  percent  difference  in  elevation -angle  error  betw'een  extreme  correc¬ 
tion  values  is  a  valid  comparison  for  this  date  or  other  dates  of  comparable 
target  altitude  and  range. 

A  secondary  comparison  included  in  Table  1  is  between  Point  Mugu  and  San 
Nicolas  Island  Raobs.  Although  not  applicable  to  a  north  cross-section  com¬ 
parison  for  Tranquillon  Peak  radar,  the  data  are  applicable  for  comparison 
of  tracking  corrections  from  Tranquillon  Peak  to  the  western  or  southern 
horizons  respectively.  One  might  accept  the  thesis  that  because  of  horizontal 
homogeneity  far  off  shore  a  San  Nicolas  Island  Raob  is  more  representative 
of  the  atmosphere  for  Tranquillon  Peak  at  low  angles  to  the  west  over  the  ocean 
than  a  balloon  sounding  over  Tranquillon  Peak  itself.  In  such  case,  one  might 
employ  the  SNI  Raob  as  the  best  representation  of  the  atmosphere.  Similarly, 
the  same  can  be  said  for  the  Point  Mugu  relationship  to  Tranquillon  Peak  for 
Southern  look-angles  (in  the  isear  shore  mixing  zone).  The  Mugu  and  Vandenberg 
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TABLE  1 


COMPARISON  OF  REFRACTION  ERRCmS 
VANDENBERG  -  NORTH 
4  October  1968 


Input 

Source 

Range  Error 

Ft 

Blev.  Error 

Mils 

Bending 

Mils 

VBG  Raob'*’ 

100.8 

3.  58 

5.60 

SNI  Raob'*' 

104.0 

2.  67 

4.70 

Mugu  Raob 

102.3 

3.  92 

6.02 

TPK  Spiral-^ 

99.4 

3. 78 

6.04 

Cross  Section 

99.6 

4.40 

7.20 

TPK  Splr/Mugu  Raob  SNI  Raob/Sect.  SNI  Raob/Sect. 
Extremes  99.4  -  102.3  2.67  -  4.40  4.70  -  7.20 

%  Difference  3  65  53 


Conditions: 

Elevation  angle  used  = 

Target  height  =  12, 000  ft. 

"■Located  directly  over  Tranq.  Peak  (TPK) 
^As  applied  to  TPK  radar 


Source 


Rauge  Error 
Ft 


Elev.  Error 
Mils 


Bending 

Mils 


r . . .  . rl 

Mugu  Raob 

Not  available  for  this  time  frame 

Mugu  Splr  1* 

86.7 

2.75 

4.  22 

Mugu  Splr  2^ 

89.8 

3.97 

5.  33 

Mugu  Section 

110.0 

3.50 

6.00 

r-—  ^  ^ 

Splr  2/Sect. 

Splr  1/Sect. 

Splr  1/Sect. 

Extremes 

86.7  -  110.0 

2.75  -  3.97 

4.22  -  6.00 

%  Difference 

27 

44 

42 

Conditions: 

Elevation  angle  =  1° 

Target  height  =  10, 500  ft 
^Located  8  am  west  of  Point  Mugj 

•4- 

Lo.?.ated  18  nm  south  of  Tranq  Peak,  75  nm  west  of  Point  Mugu 


TABLE  3 


COMi’AR ISON  OF  REFRACTION  ERRORS 

VANDENBERG  -  250®  ' 

8  June  1967 

IntMit 

Range  Error 

Elev.  Error 

Bending 

Source 

Ft 

Mils 

Mils 

VBG 

86.9 

2.97 

4.54 

Splr  2* 

89.8 

3.97 

5.33 

Splr  3"^ 

26.0 

4.31 

5.  53 

TFK  Section 

90.2 

1.  80 

3.20 

Spir  3/N.  Sect. 

W.  Sect/Spir  3 

Sect/Spir  3 

Extremes 

86. 0  -  90.  2 

1.80  -  4.31 

3.20  -  5.53 

%  Difference 

5 

140 

73 

Conditions: 

Elevation  angle  =  1°  Table  3  compares  the  pos¬ 

sibilities  open  for  Tranq  Peak 

Target  height  =  10,  500  ft  refraction  correction  for 

polar  orbit  shots. 

♦.18  nm  south  of  Tranq  Peak 
^57  nm  west  of  Tranq  Peak 


TABLE  4 


COMPARISON  OF  REFRACTION  ERRORS 
VANDENBERG  -  WEST 
30  March  1967 


.  Input 

Range  Error 

Elev.  Error 

Bending 

Source 

Ft 

Mils 

Mils 

VBG  Raob 

49.0 

1.10 

2. 17 

Splr  1* 

49.9 

1.24 

2.25 

Splr  2* 

49.9 

0.79 

1.95 

TPK  Section 

65.0 

1.60 

3. 10 

Raob/Sect. 

Splr  2/Sect. 

Splr  2/Sect. 

Extremes 

49.0  -  65,0 

0.79  -  1.60 

1.95  -  3. 10 

%  Difference 

33 

103 

59 

i 


Conditions: 

Elev.  angle  =  1° 

Target  height  -  7, 000  ft 

•Located  30  nm  and  15  ntn  west  of  TPK  coast,  respectively 
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Input 

Source 


Range  Error 
Ft 


Elev.  Error 
Mils 


Bending 

Mils 


Spir  1* 

67.7 

1,96 

! 

3.19 

Spir  3"^ 

63.0 

1.08 

2.33 

VB3  Raob 

66.2 

1.28 

2.55 

VBG  Section 

64.6 

1.60 

3.20 

Spir  3/Spir  1 

Spir  3/Splr  1 

Spir  3/Sect. 

Extremes 

59. 5  -  67.7 

1.08  -  1.96 

2.33  -  3.20 

%  Difference 

14 

82 

38 

Conditions: 

Elevation  angle  -  1  degree 
Target  height  -  8, 000  ft 
Radar  height  =  2,  159  ft 
*45  nm  south  of  Tranq  Peak 
8  nm  west  of  Point  Mugu  (within  mixing  zone) 
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TABLE  6 


COMPARISON  OF  REFRACTION  BRRCMIS 
SAN  NICOLAS  ISLAND  -  EAST-WEST 
24  July  1967 


Input  Range  Error  Elev.  Error  Bending 

Source  Ft  Mils  MUg 


SNI  Raob 

47.2 

2.58 

4. 25 

Splr  n* 

49.0 

3.68 

4. 40 

W.  Cross  Section 

51.7 

3.70 

6.60 

E.  Cross  Section 

62.4 

7.80 

9.70 

SNI  Raob/Sect.  SNI  Raob/Sea.  SNI  Raob/Sect. 

Fy-xemes  47.2  -  62.4  2.58  -  7.80  4.25  -  9.70 

%  Dlitercnce  32  202  128 


Conditions; 

Elevation  angle  =  1° 

Target  height  -  4, 200  ft  (Inversion  ton  goes  to  4, 000  ft  In  places) 
Radar  height  =  860  ft 

*Spiral  located  between  Anacapa  and  San  Nicolas  Islands. 

Ray  traps  for  these  measurements  at  6.  14  mils. 
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data  presented  show  the  corrections  are  reasonably  close.  It  Is  contrary  to 
this  reasoning  that  SNl  proHuces  less  of  a  bending  correction  than  shore,  or 
near  shore,  measurements^  such  as  displayed  in  Table  2.  The  errors  should  in¬ 
crease  to  seaward  as  the  inversion  stabilizes  and  the  marine  air  dry  air  interface 
becomes  very  sharply  defined.  This  will  be  confirmed  by  the  Table  7  sum¬ 
mary  and  will  show  also  that  full  cross  sections  to  seaward  produce  even 
greater  corrections  than  the  seaward  spirals,  except  for  the  Table  3,  June 
8,  Vandenberg  section.  Figure  16.  This  8  June  section  is  very  similar  to 
i:hat  of  29  March,  Figure  10,  with  complete  washing  out  of  the  gradient  near 
the  coast.  One  can  see  in  Figure  16.  that  the  one-degree  ray  passes  well 
above  the  sharp  inversion.*  having  penetrated  in  the  well- mixed  zone  closer 
to  the  coast,  thus  explaining  the  exception.  Recall  how  the  cloudless  zone 
of  Figure  9  and  the  hydraulic  jump  of  Figure  7  define  this  inversion  destruct¬ 
ion  phenomenon. 

In  Table  4,  Spiral  2  is  least  like  the  output  from  the  cross  section.  Spiral  2 
is  located  in  the  coastal  mixing  zone  (well  mixed  from  mechanical  eddying  ! 

induced  by  the  irnountainous  coast  itself).  Spiral  1,  30  miles  to  sea.  is  most 
like  the  cross  section. 

!  Again  hi  Table  5,  the  more  seaward  measurements  produce  the  greatest 

bending  error  v/ith  the  8 -mile  spiral  (within  the  well- mixed  zone  on  this 
occasion)  prodveing  the  smallest  errors. 

Table  6  presents  a  comparison  between  Raob,  Sprial,  and  Cross  Section  for 
SNI.  There  is  only  a  3  percent  difference  in  bending  between  the  balloon 
data  and  aircraft  spiral  data.  In  this  case,  the  spiral  is  well  to  seaward  and 
would  not  be  expected  to  differ  greatly  from  the  SNI  balloon  (Figure  11),  which 
can  be  seen  to  penetrate  the  dynamically  sharpened  inversion.  The  cross- 
section  tracing  differs  from  the  balloon  data  by  128%  due  to  the  effects  of 
steering  which  are  included  in  the  cro&s-section  correction.  Note  also  a 
considerable  diffei  ence  in  ray  traces  through  a  west  section  and  an  east  sect¬ 
ion  of  the  disturbed  atmosphere.  To  add  to  the  variation,  a  weak  split  in¬ 
version  was  found  at  many  times  (0315  L'.’,  SNI  Raob)  and  places  (Spiral  2) 
on  this  day.  The  east-west  section  of  Figure  11  shows  the  upper  inversion 
(between  ^e  290  and  260  refractopleths)  dilated  in  the  region  of  the  islands 
"bow  wave".  There  is  a  short  contraction  of  refractopleths  on  the  island's 
downslope,  followed  by  a  slight  jump  two  miles  downstream  of  the  island 
(see  Froude  number  explanations  of  flow  characteristics  in  Section  II).  The 
upper  inversion  sharply  contracts  again  downstream  of  the  hydraulic  jump. 

Here  the  one-degree  ray  receives  a  second  harsh  bending,  which  it  does  not 
receive  in  the  west  section,  thus  the  difference  in  bending  errors  in  Table  6 
between  east  and  west  sections. 

Table  7  is  a  summary  of  Tables  1-6.  Since  there  is  present  a  limited  amount 
of  statistical  items,  the  data  only  suggest  trends.  The  cross-section  data 
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TABLE  7 

SUMMARY  OF  TABLES  1  -  6 


Date 

Range  Error 
Extremes 

Blev.  Error 
Extremes 

Bending  Error 
Extremes 

4  Oct  68 
%  Difference 

Spir/VBG  Raob 

3.1 

SMI  Raob/N.  Sect. 
64.8 

SMI  Raob/N.  Sect, 
53.0 

8  June,  Mugu 

Spir  1/W.  Sect. 

26.9 

Spir  1/W.  Sect. 
44.3 

Spir  1/W.  Sect. 
42.4 

9  June,  VBG 

Spir  3/W.  Sect. 

4.9 

W.  Sect.  /Spir  3 

139.  5 

W.Sea/Splr  3 
73.0 

30  Mar  67 

VBG  RaoiyW.  Sect. 

32.7 

Spir  2/W.Sect. 
102.6 

Spir  2/W.Sect. 
59.0 

29  Mar  67 

Spir  3/Spir  1 

13.8 

Spir  3/Spir  1 

81.  5 

Spir  3/W.  Sect, 
37.8 

24  July  67 

SNIRaotj^W.Sect. 

32.0 

SNIRaoh»^W,Sect. 

202.3 

SNIRaob'W.Sect. 

128.0 

Weighted  Ave. 
%Di£f. 

17* 

107 

68 

Spiral  positions: 

4  Oct  68, 

VBG;  at  Tranq.  Peak 

8  June  67; 

8/270°  Mugu 

8  June  67 ; 

57/270°  VBG 

30  Mar  67; 

15/270°  VBG 

29  Mar  67; 

Spir  1,  45/180°  VBG 

29  Mar  67;  Spir  3,  8/270®  VBG 
■'•About  24%  when  considering  the  4  section  cases  only. 

29 

C 

produces  the  greatest  correction  requirements  for  range  error  in  4  out  of  6 
cases  by  24  percent.  The  wei^ted  average  separation  of  extremes  is  17  per¬ 
cent  The  cross  sections  produce  an  average  of  64  percent  greater  bending 
error  in  5  out  of  6  cases.  The  average  separation  of  extremes  is  68  percent. 
There  is  a  degree  of  confirmation,  though  limited,  that  Raobs  and  near  coas¬ 
tal  aircraft  spirals  (<15  nm)  produce  less  error  than  seaward  measurements. 
The  atmosphere  along  the  coast  is  better  mixed.  If  the  inversion  is  present 
it  is  split  or  dilated  and  higher  in  altitude.  As  could  be  expected,  the  greatest 
percent  differences  are  found  in  the  elevation  errors,  which  are  also  the  most 
critical  for  single  radar  target  positioning. 

Two  shortcomings  of  this  analysis  are  readily  apparent  in  this  data.  One  is 
the  small  number  of  statistics.  A  greater  number  of  cases  is  certainly  re¬ 
quired  before  the  argument  can  be  called  well  founded.  Support  for  the  mag¬ 
nitude  of  variation  is  available  in  the  work  of  Rowlandson  and  Herlihy  <1968). 
Compaxdng  aircraft  spirals  only  they  find  time  and  space  variations  of  the 
order  of  .  5  mils  for  1''  antenna  elevation  angle  when  a  series  of  spirals  are 
compared,  or  seasonal  groupings  are  compared.  The  second  a.id  more  ser¬ 
ious  shortcoming  is  that  comparisons  have  to  be  made  against  each  otner. 

An  absolute  reference  against  which  all  techniques  of  corvection  formnla 
development  and  data  gathering  systems  can  be  compared  is  urgently  needed 
in  this  type  investigation.  To  this  end  we  progress  to  the  following  bk  ?k  of 
work  on  radio  holes,  which  can  be  made  to  serve  as  this  reference.  It  will 
be  shown  that  the  cross-section  data  inputs  are  closer  to  reality  by  a  sig.iif  • 
leant  amount  (2x)  and  would  produce  the  greatest  accuracy  for  low  angle 
tracking. 
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SECTION  IV 


RADIO  HOLE  fflSSIONS 


The  nature  of  this  investigation  essentiaHy  revolves  around  the  problem  of 
bending  errors  as  they  are  affected  by  atmospheric  anomalies.  Fortimately, 
range  errors  suffer  less  deviation  than  bending  errors,  as  they  are  solely 
dependent  on  retardation  of  wave  travel.  Their  travel  time  through  the  an¬ 
omalies  is  small  compared  to  total  travel  time.  Bending  errors,  on  the 
other  hand,  respond  to  the  sharpness  of  gradient  between  the  anomalies  and 
the  ambient  condition  which  produces  a  differential  in  travel  speeds  for  por¬ 
tions  of  a  wave  traveling  through  the  gradient.  The  small  retardation  error 
means  that  a  target  can  be  positioned  within  a  very  small  triangle  of  error 
in  the  horizontal  plane,  using  ranges  from  a  radar  triangulation  system. 
However,  the  vertical  bending  of  rays  is  sufficiently  great  at  low  tracking 
angles  such  that  accuracy  of  position  in  the  vertical  plane  goes  out  of  toler¬ 
ance.  If  one  assumes  a  target  to  be  located  within  acceptable  accuracy  in 
the  horizontal  plane,  and  the  target  itself  supplies  an  accurate  measure  of 
altitude,  a  position  is  obtained  which  can  serve  as  absolute  reference  for  the 
analysis  of  bending  errors.  An  aircraft  with  C-band  beacon  using  a  radar 
altimeter  can  perform  such  a  simulation.  The  aircraft's  down-facing  radar 
altimeter  is  not  affected  by  the  anomalies  in  the  context  of  this  problem,  and 
can  be  considered  accurate  to  within  40  ft.  If  this  beacon-equipped  aircraft 
flies  to  seaward  above  the  anomalous  condition  and  experiences  Loss  of  Sig¬ 
nal  (LOS),  or  severe  signal  distortion,  allegedly  due  to  anomalous  conditions, 
then  some  combination  of  refraction  correction  formulae  and  data  inputs 
should  be  able  to  predict  the  position  of  this  signal  loss  better  than  others. 

Since  it  makes  little  difference  which  correction  formula  is  used,  concentra¬ 
tion  in  this  analysis  is  centered  on  the  role  played  by  possible  data  inputs 
and  a  consideration  of  ray  steering.  We  intend  to  show  which  of  the  already 
presented  (Section  III)  data  collection  techniques  provides  the  most  valid 
measurements  of  the  atmosphere.  If  the  measurements  have  produced  the 
right  amount  of  ray  bending,  they  should  identify  with  radio  holes  actually 
found  in  the  field;  that  is,  the  rays  should  go  negative  (bend  back  to  earth) 
precisely  where  signal  is  lost  on  the  simulation  target.  Or,  at  least  one  set 
of  measurements  should  make  a  better  prediction  of  signal  loss  than  the  others. 
A  short  explanation  of  radio  holes  and  ray  steering  is  inserted  here.  The 
time  and  position  occurrence  of  these  radio  holes  is  to  serve  as  our  absolute 
reference. 

Figure  15,  by  Ming  Wong  (1951),  shows  diagrammatically  the  causes  of  a 
radio  hole.  Observe  the  series  of  ray  traces  at  different  elevation  angles, 
from  a  transmitter  at  4,  000  ft  elevation.  Note  that  an  inversion  layer  of 
Z)  N  =  60,  is  located  at  5,  000  ft  elevation.  As  the  grazing  angle  of  the  propa¬ 
gated  wave  decreases,  there  comes  a  point  at  which  the  wave  front,  represented 
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as  rays,  perpenaicuiar  to  tne  iront  oenas  bacK  into  ine  aUuu^puex  e  uvlOW  u*c 
inversion  layer,  thus  producing  the  gap  occupying  the  center  of  the  graph. 

This  ic  due  to  retardation  of  the  lo%yer  portion  of  the  T?avc  in  moist  air  and 
velocity  increase  of  tiie  upper  portion  in  the  dryer  air.  As  a  consequence  of 
the  bending,  ducting  (or  trapping)  occurs  within  the  inversion  layer,  which 
acts  as  a  wave  guide,  and  the  signals  are  propagated  beyond  the  normal  radio 
horizon.  This  shows  up  in  Ming  Wong's  diagram  as  the  lateral  continuation 
of  waves  at  the  5,  000  ft  level.  The  trapping  can  be  advantageous  for  acquir¬ 
ing  a  target  over  the  horizon,  but  presents  a  problem  for  a  target  in  the  radio 
hole.  The  significance  of  this  loss  of  track  is  that  the  tracking  radar,  no 
matter  how  short  the  period  of  loss,  must  re -acquire  the  vehicle  in  flight. 
Recalling  that  missiles  are  traveling  on  the  order  of  thousauids  of  feet  per 
second,  decisions  regarding  missile  destruction,  or  the  attempt  at  re-acquis¬ 
ition  must  be  made  in  near  real-time.  There  is  also  a  significance  to  low 
flying  attack  aircraft  which  is  self  evident.  (Tabulated  and  diagrammatic 
compilations  of  trapping  conditions  in  general  have  been  compiled  by  Ming 
Wong,  ibid;  Stanfield,  1962:  and  Kauper,  1965. ) 

It  has  been  custom£a*y  in  the  past  to  assume  the  atmosphere  to  be  horizon¬ 
tally  homogeneous,  with  simple  stratification  in  the  vertical  plane.  The 
stratification  is  a  manifestation  of  inversion  layer  presence.  Cunningham, 
Vickers,  Fain  (1963)  and  Cunningham,  Myers  and  Crane  (1964)  attempted  to 
define  and  ray  trace  through  horizontally  anomalous  conditions  (inversions 
and  cumulus  clouds)  at  Cape  Kennedy.  The  present  work  specifically  treats 
the  Pacific  inversion.  For  this  purpose,  we  introduce  ray  steering,  as  dis¬ 
tinct  from  currently  employed  ray  tracing,  to  account  for  tilting  and  strong 
local  variation  of  gradient.  The  rays  were  steered  through  the  cross  sections. 
To  do  this,  bending  as  a  function  of  refraction  gradient  and  increasing  hori¬ 
zontal  distance  was  extracted  from  the  cross  sections  and  combined  with  the 
earth  curvature  correction  to  produce  the  corrected  elevation  at  successive 
points  in  space.  The  elementary  increments  were  chosen  sufficiently  snaall 
and  varied  to  respond  to  the  scale  sizes  of  the  anomalies  with  which  we  worked. 
A  simplifying  mechanism  and  distinctive  feature  of  the  treatment  is  the  em¬ 
ployment  of  MERS  (Milli  -  Earth  -  Radii)  as  a  fundamental  unit  of  measure¬ 
ment.  Details  are  discussed  more  fully  in  Section  IV.  The  significance  of 
trapping  is  fully  realized  in  the  case  studies  presented  below. 

.Tune  8,  1967,  Radio  Hole  Mission 


Figure  16  shows  the  cross  section  for  a  westerly  look-angle  from  Point  Mugu 
on  8  June  1967.  We  have  previously  covered  the  description  of  the  inversion, 
and  will  only  call  your  attention  to  its  upward  warping,  dilation,  lensing,  and 
consequent  weaker  gradient  toward  the  beach:  and  its  convergence  and  some¬ 
what  undulatory  nature  to  seaward.  These  are  features  that  have  strong  in¬ 
fluence  on  the  refraction  problem.  The  wave-like  convergence  at  70-80 
miles  produces  severe  bending  for  low-angle  rays. 
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It  is  this  bending  that  we  wish  to  examine.  Checking  the  originsd  radar  re¬ 
cord.  Figure  17,  the  z  plot  of  the  aircraft  climb  and  descent  shown  on  the 
cross  section  between  62  and  102  miles  is  clearly  recoi-ded  with  correlating 
time  hacks.  As  the  aircraft  descended  the  radar  elevation  angle  (z  plot) 
began  to  weaken.  At  1828  the  z  plot  of  Figure  1?  indicatpa  about  8.  200  ft 
uncorrected  for  refraction.  The  radar  altimeter  indicates  5,  400  ft  at  this 
time.  The  1830  radar  altimeter  position  of  the  aircraift  is  also  fixed  in 
Figure  16,  against  which  can  be  compared  the  1830  rays  produced  by  cross 
section  and  spiral  data.  Since  no  Raob  is  available  for  Point  Mugu,  a  spiral 
8  miles  off  the  coast  is  used,  which  should  be  a  faithful  substitution.  The 
refraction  uncorrected  ray  is  also  shown.  The  radar  at  1830  is  near  3.  4 
mils  antenna  elevation.  The  cross-section  data  indicates  the  aircraft  to  be 
3,  000  ft  higher  than  its  true  position,  and  the  spiral  (Raob)  data  indicates 
the  plane  to  be  S,  400  ft  higher.  (It  is  tme  that  these  numbers  are  not  accur¬ 
ate.  )  When  the  radar  goes  unstable,  it  is  difficult  to  pick  am  instant  in  time 
and  assign  an  elevation  angle  to  it.  This  can  be  readily  seen  in  Figure  17. 

The  Instant  of  signal  loss  represents  a  rapidly  decaying  signal  to  noise  ratio. 

The  acceptance  level  of  3  db  above  noise  is  a  statistical  criterion.  Never- 
dieleas,  the  gross  differences  Indicated  on  the  diagram  show  the  order  of 
discrepancy  involved.  This  is  evidenced  in  the  approximation  to  radio  hori¬ 
zon.  Note  the  uppermost  points  of  signal  loss  on  the  aircraft  track  slightly 
below  the  0.  0  mil  ray.  This  ray  roughly  represents  the  radio  horizon  as 
derived  through  using  cross-section  inputs  for  refraction  correction.  The 
line  of  demarcation  can  be  made  out  in  the  diagram,  except  for  the  1948  LOS 
position  which  is  very  low,  and  the  2015  AOS  which  is  very  high.  The  1830, 

2035,  1841,  2003,  and  1957  position  (from  low  center  to  upper  left  of  dia¬ 
gram)  define  a  fairly  believable  line.  The  line  suggests  that  the  cross-section 
corrections  for  the  0.  0  mil  ray  are  still  a  trifle  lean.  In  theory,  a  properly 
corrected  0.  0  mil  ray  should  lie  on  the  points  of  signal  loss. 

Of  further  interest  in  Figure  16  is  the  data  consistency  of  the  3 -mil  ray 
(trapped  in  the  inversion).  We  know  that  the  aircraft  suffers  an  increasingly 
poor  signal  starting  from  about  1828  on  the  z  plot.  Somewhere  between  1829- 
1831  the  z  and  x-y  plots  indicate  loss  of  lock.  Around  1829.  5  the  z  plot  indi¬ 
cates  about  3.  4  mils,  as  close  as  we  can  estimate.  According  to  our  ray 
trace,  3.  4  mils  looks  like  the  last  of  untrappeJ  rays.  At  3  mils  (very  close 
to  our  actual  LOS  elevation  angle)  severe  bending  takes  place.  The  3 -mil 
ray  becomes  trapped.  It  is  ducted  westward  and  just  misses  the  short  lock 
acquired  at  1917  as  the  aircraft  came  through  the  gradiert  layer.  It  is  unusual 
to  see  this  ray  steer  its  way  through  the  actually  measured  sections  of  live 
signal  between  the  100  and  120  mile  points.  The  fact  that  it  misses  the  open¬ 
ing  at  180  miles  by  only  200  meters  is  not  very  disconcerting.  The  order  in 
the  data  indicates  that  from  the  extremely  complex  set  of  measurements  a 
sensible  body  of  logic  emerges  from  which  significant  assessments  can  be  made. 

The  data  of  the  8  June  mission  suppport  the  thesis  that  Raobs,  or  near-shore 
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spirals  are  not  the  best  way  to  define  the  atmosphere  for  refraction  correc¬ 
tion.  We  conclude  in  general  that  in  the  vicinity  of  0.  0  mils,  at  100  nm,  the 
cross-section  corrected  ray  is  high  by  about  ?,  200  ft.  The  Raob  spiral  is 
high  by  4,  800  ft.  With  no  further  effort  than  the  ^iresent  project,  we  can 
state  that  low-angle  refraction  errors  can  be  halved  by  flying  cross  sections. 
Further  experimentation  could  well  improve  this. 

March  29,  1967,  Radio  Hole  Mission 

The  plan  for  the  radio  hole  missions  called  for  the  pilot  to  be  informed  by 
radio  when  a  particular  radar  lost  lock.  The  aircraft  would  continue  on  to 
seav/ard  just  above  the  inversion  until  there  was  an  acquisition  of  signal, 
that  is,  the  aircraft  would  fly  through  the  radio  hole.  When  signal  was  re¬ 
acquired,  the  aircraft  would  go  into  a  down  spiral  through  the  gradient  layer, 
and  after  passing  through  the  gradient  layer,  pick  up  a  return  course  toward 
the  beach  employing  a  porpoising  maneuver  up  through  the  layer  and  back 
down  again.  This  would  continue  until  the  radar  operators  informed  the  pilot 
that  the  plane  had  come  back  to  its  original  loss  of  signal  position  on  the 
shoreward  side  of  the  hole.  At  that  point  it  would  go  into  an  upward  spiral, 
again  sounding  the  gradient  layer  on  the  way  through.  The  upshot  cf  this  was 
that  the  aircraft  got  very  fine  detail  with  adequate  data  points  in  the  immediate 
vicinity  of  the  loss  of  lock.  What  we  hoped  to  learn  was  what  the  inversion 
layer  looked  like  (geometry  or  gradient  intensity)  at  the  immediate  point  of 
signal  loss. 

Figure  18  is  an  actual  flight  track  which  follows  this  plan.  It  is  of  interest 
to  proceed  around  the  circuit  to  show  the  correlations  br'tween  radar  signals 
and  the  structure  of  the  inversion.  First,  notice  tbe  cyclical  wave-like  course 
plotted  by  the  SNI  radar  in  the  vicinity  of  VBG  radio  hole  (2141:04  -  2144:52). 

We  have  superposed  a  copy  of  the  rcfractome.er  tra.ee  above  the  indicated 
radio  hole  on  the  diagram.  It  is  of  interest  that  the  spacing  of  waves  as  in¬ 
dicated  by  the  rofractometer  as  the  aircraft  holds  a  steady  2,  900  ft  altitude 
is  coincident  with  the  general  size  of  the  distortions  in  the  SNI  radar.  It 
follows  logically  that  while  sighting  obliauely  across  a  wave  crest  at  a  target 
crossing  the  crests  at  near  crest  height  an  azimuth  variation  can  be  produced 
between  times  when  the  line  to  the  target  is  wholly  within  the  marine  air  medium, 
and  those  times  when  it  emerges  near  the  target  in  dry  air  troughs  between  the 
crests.  This  may  explain  the  apparent  correlation  between  wave  crests,  re¬ 
corded  by  refractometer,  and  track  variation  recorded  by  SNI.  The  speculative 
quality  of  this  statement  is  emphasized. 

Proceeding  to  the  descending  spiral  left  turn  at  the  end  of  the  west-bound  leg 
(Figure  18).  the  top  of  the  gradient  layer  wa.s  penetrated  at  47:16.  SNT  lost 
track  16  seconds  before  the  gradient  layer  was  penetrated.  Note,  however, 
that  Point  Mugu  weakened  for  only  16  seconds  (4  data  fioinls).  then  came  back 
on  strong  as  the  inversion  layer  was  j'lenetra  cd,  suggesting  that  the  inversion 
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was  carrying  a  strong  Point  Mugu  sigTiai  witliin  its  "wave  guide".  Signal 
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1430,  and  did  not  give  a  reliable  position  until  the  bottom  of  the  gradient  layer 
WBB  nBRRed  through  orpoiRely  at  4R:3fi.  The  aireraft  stayed  below  the  inver¬ 
sion  with  weak,  intermittent  and  disagreeing  signals.  The  aircraft  hit  the 
bottom  of  its  descent  at  approximately  51:28,  then  began  climbing  at  that  point 
imtil  the  base  of  the  gradient  layer  was  intercepted  precisely  at  52:25. 


The  actual  instrument  traces  from  this  penetration  are  presented  in  Figure 
19.  Note  in  Figure  18  that  as  the  base  of  the  gradient  layer  was  approached, 
8dl  radars  disagreed.  As  upward  penetration  of  the  layer  progressed,  the 
tracks  merged.  Precisely  at  the  top  of  the  layer,  VBG  re-acquired,  52:52 
(compare  Figures  18  and  19  for  this  time),  but  then  drifted  into  the  same 
problem  of  short  ranging  that  it  had  on  the  opposite  side  of  the  flight  pattern 
(21:39:58  -  41:04).  The  aircraft  again  initiated  a  plunge  through  the  gradient 
and  came  through  the  base  of  the  gradient  layer  precisely  at  54:20,  when  VBG 
came  back  strongly.  The  aircraft  remained  beneath  the  gradient  layer  to  the 
end  of  the  pattern,  where  an  upward  spiral  was  performed. 


We  see  in  this  exercise  that  both  azimuth  and  range  appear  to  be  affected 
more  severely  than  previously  thought  by  the  presence  of  the  inversion  grad¬ 
ient  layer.  This  assumes,  of  course,  that  one  is  willing  to  agree  that  two 
radar  tracers  nearly  superposed,  one  on  the  other,  recording  a  strong  signal 
are  indeed  giving  an  accurate  position  in  the  lateral  plane  (with  the  exclusion 
of  the  100  ft  or  so  range  error  which  is  within  our  normally  accepted  toler¬ 
ance).  The  inter-relation  of  these  recorded  event.s  strongly  suggests  cause 
and  effect,  and  gains  a  degree  of  authority  for  lack  of  alternate  explanations. 


Let  us  now  look  at  the  atmosphere  in  cross  section  along  the  average  bearing 
from  Vandenberg  diagonally  through  the  flight  pattern  from  northeast  to  south¬ 
west.  Using  data  from  Figure  18  and  several  additional  aircraft  passes  in  this 
immediate  vicinity.  Figure  10  has  been  constructed.  Data  point  density  is 
sufficient  for  us  to  establish  that  the  strong  neck-like  convergences  of  the 
gradient  layer  and  its  general  wave  shape  are  real.  On  several  occasions 
(for  example,  the  already  presented  June  8  case)  our  cross  sections  show  that 
loss  of  lock  lias  been  coincident  with  these  wave-induced  convergences.  In 
the  example  presently  being  discussed,  the  pattern  flown  was  the  second  of 
two  flown  at  the  exact  same  point  for  the  same  radio  hole,  spread  in  time  by 
80  minutes.  It  appears  that  these  anomalies  can  remain  as  standing  wave- 
like  phenomena  for  considerable  time.  Recall  also  the  time  variation  pre¬ 
sented  in  Figure  11  of  the  shift  of  wave  nodes  over  Tranquillon  Peak  during  a 
45-minute  interval.  Apparently,  at  sea  the  waves  are  considerably  more 
stationary,  sharply  defined,  and  of  intense  gradient;  as  opposed  to  the  more 
diffuse  topographically  induced  waves  found  over  Tranquillon  Peak. 


Consider  the  cross-section  revelation  of  bending  effects.  We  have  ray-traced 


through  the  cross  section  of  Figure  10  using  the  three  possible  data  inputs 
which  we  have  previously  discussed,  raobs,  spirals,  and  cross  sections. 

In  this  instance  we  have  more  precision  than  that  obtained  for  our  previous 
8  June  example.  For  this  March  29  instance,  there  is  excellent  agreement 
between  the  aircraft  position  for  loss  of  lock  (41:04,  Figure  3  8)  and  the  an¬ 
tenna  elevation  at  the  time  of  signal  loss.  That  is,  agreement  is  good  when 
the  antenna  elevation  is  corrected  for  refraction  by  cross-section  inputs. 

Rays  for  balloon  corrected  and  aircraft  spiral  corrected  aitenna  elevations 
are  shovm  for  comparison.  Again,  the  improvement  of  cross-section  data 
inputs  coupled  with  ray  steering  calculations  defined  the  most  accurate  path 
to  low -angle  targets. 

Mission  Abort,  or  Missile  Destruction  Condition 

According  to  Range  Operations  personnel,  signal  strength  must  be  3  db 
above  noise  for  silo  missiles  to  check-out  before  launch.  Failure  to  obtain 
this  minimum  ratio  is  cause  for  mission  abort.  Moreover,  failure  of  radars 
to  agree  on  missile  position  immediately  after  launch  is  cause  for  destruction 
of  the  missile.  Let  us  see  how  such  negative  situations  can  come  about  as 
direct  effects  of  atmospheric  anomalies.  Recall  Figure  7,  which  shows  how 
th'.'  inversion  layer  can  present  severe  trapping  hazards  for  Tranquillon  Peak 
radar.  We  have  ray-traced  through  this  section  and  find  that  the  rays  will 
trap  below  50  mils,  given  these  conditions.  Moreover,  it  is  possible  for  the 
main  lobe  of  the  propagation  pattern  to  be  trapped  while  the  side  lobe  locks 
onto  a  missile  coming  up  through  the  inversion  from  its  silo.  Assuming  0.  8° 
minimum  lobe  separation,  the  altitude  error  is  on  the  order  of  1500  ft  at  18 
miles  relative  to  the  main  lobe.  For  a  distant  radar  operating  under  similar 
conditions,  such  as  Point  Pillar,  tracking  a  Vandenberg  Minuteman,  the  al¬ 
titude  disagreement  between  VBG  and  Point  Pillar  radars  would  be  greatly 
increased.  Point  Pillar  could  experience  inversion  "lobe  splitting"  just  as 
easily  as  Tranquillon  Peak,  with  disagreement  against  Tranquillon  Peak's 
missile  altitude  running  around  6,000  ft.  We  have,  therefore,  all  the  in¬ 
gredients  needed  to  produce  either  a  mission  abort  due  to  weak  signal,  or  no 
signal;  or  a  mission  destruct  due  to  confused  radar  after  the  missile  has  left 
the  pad. 

Summary 

We  have  presented  three  radio-hole  situations  which  have  provided  either 
direct  confirmation  or  logical  explanations  of  the  following  effect: 

1)  full  cross  sections  employing  steering  corrections  will  produce 
greater  accuracy  for  refraction  correlations  by  a  factor  of  2,  than  correc¬ 
tions  based  on  raob  or  aircraft  spiral  data  sources; 

2)  azimuth  variation  appears  to  be  more  strongly  affected  by  stratifi- 
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cation  irregularities  than  previously  believed; 

3)  severe  bending)  trapping)  and  lube  splilUiig  due  Lo  vai'lalions  uf 
the  inversion  offer  strong  logic  for  explaining  beacon  fade,  loss  of  signal, 
confused  radar  and  consequent  missile  destruction. 
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CHAPTEE  V 


EXPOSITION  ON  mathematics  EIvIFLOYeD 
IN  REFRACTION  ANOMALY  INVESTIGATION 


Sections  I-IV  have  presented  the  results  of  field  measurement  analysis.  In 
conjimction  vdth  the  field  program  significant  strides  were  made  in  the  refine¬ 
ment  of  computational  techniques.  Starting  with  a  review  of  RBFRAC  program 
(refraction  correction,  supplied  by  Western  Test  Range  Eiigineering  Depart¬ 
ment).  consideration  was  gfven  to  simplification  of  correction  techniques. 

This  was  dictated  by  needless  complexities  revealed  in  the  review  and  by  our 
own  lack  of  an  adequate  computer  to  handle  the  storage. 

Secondly,  it  was  quickly  learned  that  the  simplification  efforts  led  to  the  cap¬ 
ability  of  real-time  calculations  and  trapping  warnings  for  low  and  negative 
angles. 

Thirdly,  one  could  no  longer  ansume  the  atmosphere  to  possess  horizontal 
homogeneity  when  dealing  with  low-angle  tracking;  nor  could  atmospheric 
variations  be  ignored  on  grounds  that  the  effects  are  small;  nor  can  we  con¬ 
veniently  claim  the  required  measurements  are  impossible  to  get.  The  real 
atmosphere  demands  that  the  irregularities  be  measured  and  accounted  for. 
That  is,  signals  do  not  systematically  penetrate  through  a  well  organized 
medium,  but  must  be  steered  through  locally  tilted  surfaces,  converging  re¬ 
fraction  gradients,  and  lensed,  sandwiched,  or  washed  out  inversion  struc¬ 
tures.  To  serve  this  problem,  a  simple  ray  steering  technique  was  developed, 
and  this  program  along  with  considerations  of  the  above  Issues  is  presented 
herein. 

The  aims  of  the  present  section  are,  therefore: 

1)  Review  the  REFRAC  program  supplied  by  the  Western  Test  Range, 
and  modify  for  present  use. 

2)  Present  real-time  refraction  correction  and  trapping  warning  tech¬ 
niques  good  for  all  angles. 

3)  Acknowledge  the  refraction  complexity  of  the  atmosphere  and  cope 
with  its  irregularities  by  ray  steering  techniques. 

Review  and  Modification  of  REFRAC  Program 

The  REFRAC  program  was  reviewed  and  found  mathematically  correct,  but 
overly  complex.  We  were  unable  to  reproduce  trial  runs.  Modification  was 
required  to  cure  such  problems  as  computational  instability  and  to  meet  our 
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LIST  OF  SYMBOLS 


A  :  a  variable  in  refraetiun  prugrain 

A.  (DEN/ 2):  a  constant  in  real-time  program 

AUX  :  an  auxiliary  variable  in  refraction  program 

B  :  a  variable  in  refraction  program,  also  a  constant  in  real-time 
program 

B  :  initial  radar  elevation  angle,  initial  emergence  angle,  apparent 
target  elevation  angle 

Bj  ;  emergence  angle  at  ith  level,  angle  of  ray  with  the  horizon 

tio  :  true  elevation  of  target  as  seen  from  radar 

c  :  speed  of  light  in  a  refractive  medium 

c^  :  speed  of  light  in  vacuo 

C  :  a  variable  in  refraction  program,  also  a  constant  in  real-time 
program 

Curv,  l/Rj,:  instantaneous  curvature  of  ray 

D  :  a  variable  in  refraction  program,  also  a  constant  in  real-time 
program 

DEN  :  normalized  decay  parameter  in  dimensionless  units  in  single  ex¬ 
ponential  atmosphere  model 

DI  :  an  initial  value  in  refraction  program 

Dj,  D^:  normalized  decay  parameters  in  dimensionless  units  of  quadratic 
exponeiitial  model  atmosphere 

^li’  ^2i;  refraction  parameters  used  in  Gardner's  (1964)  REFRAC  program 
:  prefix  indicating  increment 
Ah^  :  change  in  height  of  target  at  ith  level 
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:  N-unit  change  over  fixed  interval  of  .  01  mers  in  ray-steering 
program 

^  :  bending  increment  at  ith  level 

2  2 

E,  exp(t  -t  ):  a  vEiriable  in  refraction  program,  also  a  computationsd 
parameter  in  real-time  program 

c  :  elevation  error 

y  :  auxiliary  angle  between  and  j3^  used  in  Rowlandson’s  formulation 
j  (5  •  doppler  radar  error  angle 

I  t  ^  :  first  iterated  value  of  elevation  error 

F  :  a  variable  in  refraction  program 

F(x)  :  a  special  refraction  function  in  real-time  program 
FI  :  an  initial  value  in  refraction  program 

G  :  a  variable  in  refraction  program,  also  a  computational  parameter 
in  real-time  program 

h  :  elevation  above  sea  surface 

h.  :  elevation  of  ith  level  in  stratified  atmosphere,  also  elevation  of  ith 

point  in  ray-steering  program 

h  :  radar  elevation 

0 

h^^  :  height  of  last  point  in  ray-steering  program 

H,  h/r  ;  elevation  in  normalized  dimensionless  units 
o 

HD^  :  range  in  mers  of  ith  iteration  point  in  ray-steering  program 
HD^  :  range  in  mers  of  last  point  in  ray-steering  program 
i  ■  subindex  for  ith  point 

KI,  1-cos  an  initial  value  used  in  refraction  computations 
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KB,  1-cos  B:  a  variable  used  to  compute  emergence  angle 

mer  :  distance  subtended  at  earth's  surface  by  a  central  angle  of  1  mil 

n,  c/c  ;  index  of  refraction,  also  subindex  for  nth  point 
o 

:  index  of  refraction  at  the  ith  level 

••6 

N,  10  .  (n-1):  practical  unit  of  index  of  refraction 

:  index  of  refraction,'  in  practical  units  at  ith  level 

N  :  index  of  refraction,  in  practical  units,  at  radar  elevation,  also  a 
^  refraction  parameter  in  simple  exponential  atmosphere  model 

N  :  average  index  of  refraction,  in  practical  units,  over  a  layer 

o  :  subindex  indicating  initial  conditions 

r  :  radius  vector  from  center  of  the  earth  ^ 

r.  :  radius  vector  to  ith  level 

1 

r  :  radius  vector  to  radap:  site 

o 

R  :  earth's  radius 

Rc,  1/Curv;  local  radius  of  curvature  of  ray 

R^j  ;  mean  radius  of  curvature  of  ray  in  the  ith  layer 

RGBRR:  range  error 

SE,  sin  (THE):  sine  of  angle  subtended  by  radar  range  as  seen  from  the 
center  of  the  earth 

SR  ;  radar  electromagnetic  range 

S,  SR/r^:  radar  electromagnetic  range,  normalized  to  dimensionless  units 
SR^  :  radar  electromagnetic  range  of  ith  point 

SEGA,  sec  (B^  j):  secant  of  lower  emergence  angle  of  a  computational  layer 
SECI.  sec  (B^):  secant  of  initial  radar  elevation  angle 
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SECB,  sec  (Bj):  secant  of  upper  emergency  angle  of  a  computational  layer 

TANA,  tan  (Bj_j):  tangent  of  lower  emergence  Euigle  of  a  computational  layer 
TANI,  tan  (B^):  tangent  of  initial  radar  elevation  angle 

TANB,  tsui  (Bj):  tangent  of  upper  emergency  angle  of  a  computational  layer 

TAU  :  total  bending 

THE,  angle  subtended  by  radar  range  as  seen  from  the  center  of  the  earth 
TRGE:  true  range  of  target 

t^  :  an  auxiliary  constant  in  read-time  refraction  program 
t  :  an  auxiliau*y  variable  in  real-time  refraction  program 

T  :  path  bending  amgle 

V  :  a  dummy  variable 

V  ;  shorthand  for  1  -  cos  <t) 

W,  1-cos  (THE)'  am  auxiliary  variable  in  refraction  program 
X  :  an  auxiliary  computational  variable 

0*,  THE:  amgle  subtended  by  radao*  range  as  seen  from  the  center  of  the  earth 

z  :  a  computational  variable 

2 

zz,  1/(1  -  z  ):  an  auxiliaury  variable 
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Own  restraints  of  using  a  simple  computer. 

Review  of  Pertinent  Refraction  Correction  Principles 

Ray-tracing  techniques  are  well  developed  and  we  need  not  repeat  here  any 
basic  derivations.  Essentially,  all  computing  methods  are  based  on  the  fact 
that  a  ray  traveling  in  a  refractive  medium  will  bend  so  that  its  local  curv¬ 
ature  (defined  as  the  reciprocal  of  the  instantaneous  radius  of  curvature  R  ) 
is  given  by  ^ 

Curv  =  l/R  =  1/n  (dn/dR  )  (1) 

c  c 

in  which  "n"  is  the  local  index  of  refraction.  This  is  actually  Snell's  law  in 
differential  form.  Differential  equation  (1)  can  be  solved  either  formally  or 
numerically  as  soon  as  a  distribution  for  "n"  is  specified,  since  the  radius 
of  curvature  is  by  definition  normal  to  the  trajectory. 

The  usual  approach  to  the  problem  is  to  divide  the  atmosphere  into  a  number 
of  sufficiently  thin  slices  so  that  the  vertical  gradient  of  "n"  in  each  slice  is 
sensibly  constant.  When  this  is  done,  equation  (1)  can  be  integrated  formally 
if  several  further  simplifying  assumptions  are  made  (Schulkin,  1952),  (Weis- 
brod  and  Anderson,  1959),  (Gardner,  1964). 

Given  a  stratified  atmosphere,  it  cam  be  shown  that  the  change  in  direction 
of  the  emergent  ray  between  any  two  given  elevations  depends  only  on  the 
initial  and  final  indices  of  refraction  (Bean  and  Dutton,  1966).  This  is  ex¬ 
pressed  by  the  equation 


n.r.  Cos  B.  =  n  .  r  .  cos  B  =  constant  (2) 

1X10  0  o 


in  which  B.  is  the  emergence  angle  (the  angle  between  the  ray  and  the  horizon) 
and  r,  is  t^e  vector  polar  coordinate  with  origin  at  the  center  of  the  earth. 
EquaUon  (2)  is  the  finite  form  of  SneR's  law  for  an  isotropic  spherically  strat¬ 
ified  medium. 


All  solutions  for  the  ray-tracing  problem  which  are  based  on  these  premises 
differ  only  in  the  mathematical  simplifications  invoked.  For  instance, 
Schulkin' s  (ibid)  solution  for  the  local  elementary  bending  of  the  ray  for  low 
angles  of  elevation  between  two  slices  is: 

At.  =  2  .  (N.  ,  -  N;)  /  (B.  +  B.  ,)  .  10‘®  (3) 

1  1-1  1  1  1-1 

where  is  the  refractivity  in  "N"  units, 

N.  =  (n.  -  1)  .  10®  (4) 
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v/hile  Weisbrod  and  Anderson  (ibid)  give  theirs 
At^  =  2  .  (N._j  -  N^)  /  (tan  + 

If  sufficient  slices  are  taken,  and  trigonometric  subroutines  are  used.  (3) 
will  be  faster  in  the  computer  than  (5),  while  maintaining  an  accuracy  com¬ 
patible  with  the  usually  available  inputs.  Gardner's  (ibid)  solution,  which 
uses  the  subtended  earth's  angle  0*  as  the  main  parameter,  turns  out  to  be 
computationally  equivalent  to  the  Weisbrod  and  Anderson  solution,  as  can  be 
seen  by  the  fact  that  Gardner's  modified  equation  for  the  elementary  traject¬ 
ory  length  (Gardner,  1964,  page  16) 


as 


tan  B.  ,)  .  10 
1-1 


(5) 


10 


-6 


SR.  -  SR.  , 
1  1-1 


N..  (H.  -  H.  ,) 
1  1  1-1 


(sec  B.  +  sec  B.  , . 
1  1-1) 


c  .  (tan  B.  +  tan  B.  , ) 
1  1-1 


(6) 


is  equivalent  to  Weisbrod  and  Anderson's  formulation 

SR.  -  SR.  =  Rc..  At.  (7) 

1  1-1  i  1 

which  solves  for  the  elementary  arc  length  in  terms  of  the  local  mean  radius 
of  curvature,  Rc.  and  the  local  bending,  since  the  mean  radius  of  curvature 
Rc.  is  ^ 

5c.  10°  ,,, 

^  c  .  2  .  (N.  -  N.  J 

1  1-1 

Introducing  (8)  and  (5)  into  (7)  actually  reproduces  (6). 

Gardner's  solution  of  the  ray-tracing  problem  when  a  constant  gradient  is 
assumed  for  the  refractive  index  between  layers  is  unnecessarily  complex 
computationally.  For  example,  differences  between  small  and  very  similair 
quantities  are  often  programmed,  as  in  his  proposed  equations  (29)  and  (31), 
(Gardner,  l'J64) 


SR.  -  3R.  , 
1  1-1 


c  . 

■  __P 
c 

tan  B 


D„.  .  cos  B  .  (tan  B.  -  tan  B.  ,) 
2i  o  1  1-1 

cos  0  -  -  r  /  r" 

_ _ o 

sin  6''- 


which  although  quite  correct  mathematically,  can  lead  to  computational  in¬ 
stability  at  low  angles  of  elevation  due  to  cancellation  of  leading  digits. 

When  using  a  high-speed  computer  to  solve  a  problem  it  is  generally  best  to 
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putaiiunal  stc-fis;  thr  computer  takes  far  Ififis^er  in  makirsfi  ingjt-ai  ili-LiaiOtin 
or  in  performing  trigonometric  or  similar  subroutines  ttian  in  (Joing  addilion-; 
or  muuipUcations.  We  iiave,  Ihereforr,  jtrivrri  to  simplify  the  mathematical 
routines  as  much  as  possible  and  to  avoid  unnecessary  storage,  always  check 
ing  the  results  against  such  standards  as  the  precomputed  tables  for  the  CRI’L 
reference  atmospheres  or  Gardner's  (ibid)  computations  to  insure  that  con¬ 
sistency  Is  maintained. 

Modification  of  Ray  Tracing  for  Simple  Computers 

Available  computing  methods  were  reviewed  and  modified  to  devise  a  simpler 
routine  which  would  insure  computational  stability  and  minimum  storage  re¬ 
quirements.  The  programming  was  arranged  to  avoid  as  many  trigonometric 
subroutines  as  possible.  In  the  end  we  selected  Weisbrod  and  Anderson's 
formulation  (5)  for  the  bending  rather  than  Schulkin's  (3),  as  it  is  possible  to 
program  the  former  without  using  many  trigonometric  subroutines  so  that  it 
computes  about  as  fast  as  Schulkin's  while  retaining  the  ability  to  work  with 
high  angles.  Most  trigonometric  subroutines  aind  large  storage  requirements 
can  be  avoided  in  the  follov/lng  manner. 

The  initial  value  KI  =  (1  -  cos  B  )  is  computed  accurately  using  a  rational 
approximation  designed  specifically  for  it.  Using  the  built-in  computer  tri¬ 
gonometric  subroutine  for  the  cosine  function  will  lead  to  cancellation  of 
leading  digits  at  low  angles  of  elevation 

z  =  B  .  (0.  0837984973-4.  15747751/(B^  -9.  98909688)  ) 
o  o 

zz  =  1/ (1+z^)  (9) 

KI  =  1  -  cos  B  =  2  .  z^  .  zz 
o 

The  following  initial  constants  are  then  stored  for  the  starting  level 


SECI  =  l/zz/(l  -  z  ) 
TANI  =  2  .  z/(l  -  zh 
A  =  10"?N 

o 

C  =  h  /H 
o 

F=A+C+A. C 
FI  =  F 

TANA  =  TANI 
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n 

At  the  ith  level  (i  =  1.  2.  ..  n),  the  following  quantities  are  coniputed  Ir.  succ- 
esaion  to  calculate  the  bending  and  the  electromagnetic  range 

B  -  lO'f  N. 

1 

D  =  h./R 

E  =  B+  D+  B.  D 
G  =  (E  -  FI)/(1  +  E) 

KB  -  KI  -  G  .  (  1  -  KI)  (10) 

SECB  =  1/(1  -  KB)  .  , 

TANB  =  KB  .  (2  -  KB)  )  '.  SECB 

TAU  =  TAU  +  2  .  (A  -  B)  /  (TANA  +  TANB) 

SR  =  SR  +  B  .  R  .  (D  -  C)  .  (SEC A  +  SECB) /(TANA  +  TANB) 


TANA  =  TANB 
SECA  =  SECB 

A  =  B  (10a) 

C  =  D 
F  =  B 

The  total  accumulated  bending  is  TAU,  auid  the  last  five  statements  are  a 
reinitialization  of  conditions.  We  can  now  move  to  the  next  level  and  never 
have  more  than  just  a  few  constants  in  storage.  The  only  subroutine  needed 
is  the  square  root.  Since  each  successive  KB  is  computed  from  the  starting 
value  Kl,  there  is  no  build-up  of  rounding-off  errors. 

Following  Gardner  we  now  move  to  the  (i  +  l)th  layer  and  iterate  these  cal¬ 
culations  until  the  calculated  electromagnetic  path  exceeds  the  recorded  one. 
The  final  ith  level  is  then  adjusted  by  a  proportional  interpolation  (Gardner, 
1964),  so  that  the  computed  range  coincides  with  the  observed  one.  All  the 
other  parameters  needed  for  the  ray  tracing  can  be  derived  from  the  total 
bending.  We  have  used  the  subtended  angle  at  the  center  of  the  earth  "THE" 
as  the  parameter  for  the  true  range  calculations,  but  obtained  it  from  the 
bending  rather  than  vice-versa  as  Gardner  does.  This  angle  is  obtained  from 
the  relation 

I’HE  =  (B.  -  B  )  +  TAU  (11) 

1  o 


51 


J 


i  ^  '  I  (  J  •  i  .i!  .i<  ‘  i '  -  ’  'I  :p'!'  1  'U  r  I  •  I ;  ,  !!  >  f  •  < f  •  i  .  !ii  .  •  ’  i  i 
ti.  L'V  inij*  t  If  ! I  f ff iji lij i uiM.1  i'_  f  u.> ubr ■  ifi ' '  1  i’‘, 

i  ANR  -  TAM 

-  ■  /«  .  H  :  (  1  J  1 

- . 

Siruje  al.  all  practical  radar  ranges  the  angle  (B.  -  B  )  is  small,  wc  ean  then 
compute  (he  inverse  function  arctan  (x)  hy  a  local  rsSional  approximation  de¬ 
signed  specifically  for  small  angles,  which  is  faster  than  using  the  computer 
built-in  inverse  trigonometric  function  subroutine, 

X  .  (6.  36918871/(x^+8,  60141922-4.  43698658/ (x^+ 1.  9876921)))  (12a) 

The  true  range  is  obtained  next  by  solving  the  triangle  with  vertices  at  the 
center  of  the  earth,  at  the  radar  site  and  at  the  intersection  of  the  adjusted 
idi  level  with  the  ray.  Since  the  angle  subtended  by  the  total  trajectory,  as 
viewed  from  the  center  of  the  earth,  is  small  for  all  practical  radar  ranges, 
we  found  it  preferable  to  obtain  the  value  W  =  1  -  cos  (THE)  using 

W  =  THE^.  (-0.  745585589  +  37.  235360/(29.  8938597  +  THE^))  (13) 

and  then  obtain  the  true  range  from 

TRGE  =  R  .  ((D-DI)^  +2.  D.  DI.  W)^^^  (14) 

which  avoids  computational  instability  due  to  loss  of  leading  digits.  The 
range  error  is  then  given  by  the  difference 

RGERR  =  TRGE  -  SR  (15) 

where  we  have  followed  the  usual  mathematical  convention  to  define  the  error 
as  the  difference  between  the  true  value  and  the  approximation  rather  than 
vice-versa.  An  error  is  the  amount  to  be  added  (algebraically)  to  an  obser¬ 
vation  to  bring  it  up  to  its  true  value. 

Finally,  the  elevation  error  "c  "  is  computed  as  follows 

Z  =  T.  (0.  0837984973-4.  15747751 /(t“-998909688)) 

V  -  2  .  z/(l  +  z2) 

AUX  -  (D  -  DI)/(1  +  D)  (16) 

SE  =  2  ,  z/(l  +  z2) 

X  =  (V-SE.  TANB  -  AUX)/(SE  +  TANI.  (1  - AUX)-TANB.  (1-V)) 

C  =  X  .  (0.  443525198  +  9.  28995955/(x2  +  16.  6943041)) 

The  various  approximations  used  are  valid  to  at  least  7  decimals  for  B 
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The  program  waa  completed  with  the  addition  of  safeguards  to  warn  of  poten¬ 
tial  or  actual  trapping  layers.  In  the  latter  case,  the  program  is  capable  of 
following  the  trapped  ray  back  to  earth  smd,  assuming  a  specular  reflection, 
to  follow  the  ray  for  as  many  reflections  as  desired.  If  one  or  more  poten¬ 
tial  trapping  layers  are  implicit  in  the  inputs,  the  largest  angle  of  penetra¬ 
tion  is  computed  and  printed  at  the  end  of  each  run.  The  program  accepts 
positive,  negative  or  zero  elevation  angles  of  any  magnitude.  If  the  input  is 
a  negative  elevation  angle,  the  antenna  elevation  angle  and  range  of  the  radio 
horizon  are  also  automatically  computed  and  printed.  If  the  input  elevation 
angle  is  smaller  than  the  radio  horizon  antenna  elevation,  a  specular  reflect¬ 
ion  is  assumed  where  the  ray  hits  the  earth  and  the  bounced  ray  is  followed 
through.  A  modification  of  the  computations  is  necessary  in  this  case  since 
formula  (11)  no  longer  holds.  Since  a  target  can  be  reached  by  either  a 
direct  or  a  close-by  ground -reflc -'ted  ray  at  low  elevations,  a  fact  which 
causes  no  small  amount  of  trouble  in  the  hiuiting  mechanism  of  the  radar  at 
low  elevations,  the  ability  to  analyze  both  cases  is  significant.  In  spite  of 
the  relative  progran  ''ohistication,  only  41  computing  statements  involving 
about  77  multiplicati<  j  and  divisions,  60  additions  or  subtractions,  2  square 
root  routines,  and  11  logical  decisions  per  layer  are  needed  for  the  non- 
trapped  ray.  Although  the  present  program,  written  in  TELCOMP  language, 
serves  well  the  purposes  of  the  present  investigation,  a  modest  additionsil 
effort  produces  real-ume  corrections  with  an  accuracy  compatible  with  the 
assumptions  inherent  In  stratification. 

Real-Time  Corrections  for  Low-Angles  of  Elevation 

We  have  seen  that  present  reduction  methods  share  in  common  the  assump¬ 
tion  of  a  stratified  atmosphere  in  which  the  index  of  refraction  varies  linearly 
with  height  within  selected  layers.  When  the  layers  are  not  too  far  apart, 
this  is  a  fair  assumption,  but  the  need  for  many  reference  layers,  which  re¬ 
quire  separate  individual  computations,  greatly  increases  the  execution  time. 
One  alternative  is  to  precompute  and  to  store  the  results  of  several  key  pro¬ 
files  and  to  interpolate  for  the  corrections,  but  this  requires  a  large  storage 
capacity. 

Refractivity  profiles  generally  decay  exponentially  with  height  when  one  looks 
over  fairly  thick  layers,  so  it  is  obvious  that  the  thickness  of  the  reference 
layers  could  be  increased  considerably  if  suitable  decay  constants  were  avail¬ 
able.  Freeman  (1964),  following  the  lead  of  Thayer  (1961),  has  shown  how 
acceptable  approximations  can  be  made  to  yield  a  closed-form  solution  for  the 
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Since  GGrrcctiGns  srs  greatest  at  the  '*«»  pnoles  of  elevation  which  are  of 
primary  interest  in  this  report,  and  since  the  approximations  Invoked  by 
Freeman  are  critical  in  this  area,  we  have  revised  Freeman's  derivation 
keeping  r.econd-order  terms  in  order  to  investigate  the  (lossitaiiixy  uf  im¬ 
proving  his  solution. 

Correction  Derived  From  Using  the  Curved  Radar  Path  Instead  of  the 
Geometric  Slant  Range 


Consider  a  bent  path  between  two  points  A  and  B  (Fig.  20)  for  which  the 
total  bending  is  T.  The  radius  of  curvature  changes  continuously  from  A  to 
B  but  we  can  obtain  a  good  estimate  of  the  error  Involved  in  substituting  the 
curved  path  ^8  by  the  chord  AB,  by  assuming  an  average  radius  of  curvature 
R  between  A  and  B  will  cut  the  segment  AB  at  an  angle  t/2.  The  relative 
error  committed  in  substituting  this  arc  by  the  chord  is 


_ arc  -  chord.  ,  sin  (T/2) 

Rei  err  - - - - —  -  i  -  .  i^. - 

arc  (T  /  2) 


T^/24 


(1  -  T‘‘/80 . ) 


(17) 


The  angle  T  is  generally  less  than  20  mils  or  so,  so  the  error  introduced 
by  this  substitution  is  generally  less  than  one  fifty  thousandth  of  the  range, 
A  total  bending  of  nearly  25  degrees  or  430  mils  would  be  needed  to  intro¬ 
duce  a  one  percent  range  error.  Bendings  of  this  magnitude  occur  only 
under  circumstances  approaching  the  limits  of  the  assumptions  invoked  for 
rny-traoing  We  should,  of  course,  add  this  correction,  especially  when 
dealing  with  long  ranges  and  bendings. 

Target  Elevation 


In  what  follows  we  will  consider  all  distances  normalized  to  dimensionless 
units  by  expressing  them  in  terms  of  the  length  of  the  radius  vector  r^  from 
the  center  of  the  earth  to  the  radar  site.  This  makes  it  easier  to  keep  track 
of  the  order  of  magnitude  of  the  various  terms  involved  and  simplifies  com¬ 
putations.  The  radius  vector  to  the  radar  site  is 


r  =R  +  Ah  (18) 

o  o 

where  R  is  the  radius  of  the  earth  (in  arbitrary  units)  and  is  the  elevation 
of  the  radar  site  above  the  earth's  surface  in  same  units.  The  normalized 
radar  range  S  is  then 


S  =  SR/r 

o 


(19) 
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H  Ah/r  (20) 

o 

If  jS  is  the  true  target  elevation  angle  along  the  geometric  path  SR.  and  R 
the  actual  radar  elevation  angle,  v.'e  can  obtain  from  triangle  OSK  in  Figure 
21  the  normalized  target  elevation 

H  =  (1  +  +  2  .  S  .  sin  /3  -  1  (21) 

o 

Bending  Error 

An  exponential  profile  is  next  assumed,  with  a  reference  surface  refractive 
index  N  (in  "N"  units)  and  a  decay  constant  DEN  so  that  the  refractive  index 
at  the  normalized  height  H  above  the  reference  surface  is  given  by 

N  =  N  .exp  (-DEN  .  H)  (22) 

o 

We  must  note  that  the  decay  constant  must  also  be  expressed  in  terms  of  the 

radius  vector  r  . 

o 

If  the  atmosphere  is  stratified,  the  ray  will  obey  Snell's  law  in  spherical 
coordinates  (2)  and  in  normalized  units  the  emergence  angle  B  at  elevation 
will  satisfy  the  relation 

cos  B  =  n  .  cos  B  /(n  .  (1  +  H))  (23) 

o  o 


in  which 


n  =  1  + lO”®  .  N 
o  o 

n  =  I  + lO"®  .  N 

substituting  (22)  and  (24)  into  (23)  yields 


(24) 


(1  +  10"®  .  N  ).  cos  B 
„  o  o 

cos  B  = 


(25) 


(1  +  10' 


We  will  next  assume  that  the  normalized  range  is  of  the  order  10  ^  or 
smaller.  Expanding  (21)  to  a  second  order  approximation 

H  =  1/2  .  [  (S  .  cos  +  tan  0^)^  -  tan^  (26) 

-  6 

Furthci-more.  the  quantity  10  .  N  .  exp  (-DEN  .  H)  in  the  denominator 

of  (25)  never  exceeds  5  .  10'^,  so  i? is  always  three  orders  of  magnitude 
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o  o 

which  indif-atcLs  that  far  amall  an(/les  of  clevatlori.  It  is  ai  least  one  or'ier  of 
magnitude  smaller  than  S  within  the  ansuniptions  made.  We  can,  therefore, 
to  at  lea.si  a  second  order  of  approximate cn,  rewrite  (25) 

~  (1  +  10“®  .  N  .  (1  -  exp(-UEN  ,  H)  -  II)).  cos  B  (28) 

COS  I5  -  O  O 

The  elementary  bending  at  any  point  of  the  path  is  given  by 

dT  =  (i/n)  .  (dn/dH)  .  cos  B  .  dS  (29) 

Taking  derivatives  in  (22)  and  combining  with  (24)  and  (28)  we  get 

dT=-10‘®,  N  .  DEN.  cosB  N  ,  (1-2.  exp(-DEN.  H))-H].  exp(-DEN.  H).  dS 

o  o  o 

(30) 


We  will  now  introduce  a  new  variable  to  simplify  the  handling  of  the  various 
terms.  Let 


A  =  (DEN/2)^^^ 

t  =  A  tan.  5  (31) 

o  o 

t  =  t  +  A .  S  .  cos  B 
o  o 

Then;  (26)  can  be  written 

~  2  ?  2 

H  =  (t  -  L)/2A  (32) 

and  hence 

dS  -  dt/(A  ,  cos  B  )  (33) 

o 

so  equation  (30)  becomes 

dT  ~-2.  10‘®  .N^.  cosB^.  sec^^.  (I+IO'®.  N^.  (1-2.  exp(t^-t^)+(t^^-t^)/2A^].  exp(t^-t^).  dt 

(34) 


which  can  be  integrated  formally  to  yield 

T  =  B.  (F(t  )  -  E.  F(t))+D.  (t  -E.  t)  +C.  (F(/2t  )-E^.  F/2t))  (35) 

o  o  o 
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exp 


(x**)  .  f  ^  esp(-v‘  i  .  dv 


(dii) 


The  function  F  (x)  varies  smoothly  from /ff  to  zero  when  the  argument  x 
varies  from  zero  to  an  infinitely  large  value,  and  admits  a  simple  and  accur¬ 
ate  rational  approximation.  When  second  order  terms  are  retained  in  (3jl 
the  constants  are 


B  =  10"®  .  N  .  A  .  (1  +  lO"®  .  N  +  (2t^  -  1)/4A^) 
o  o  o 

-19  9 

C  =  2  .  10  .  N  .  A 

o 

D  =  .  10’®  .  N  /2A 
o 

E  =  exp  (t^  -  t^) 


F(x)=l.  00129/(x+.  039262+.  1856633/ (x-9.  024780+115.  2736/(x+12.  29742))) 

The  rational  approximation  developed  for  F(x)  is  accurate  to  about  five 
decimals  for  all  ranges  of  x.  All  of  these  expressions  can  be  computed 
very  fast  even  on  a  small  computer. 


An  analysis  v/’as  then  made  of  the  extreme  possible  values  of  all  the  con¬ 
stants  in  <37).  The  two  correcting  terms  which  multiply  constants  C  and  D 
turned  out  to  be  at  least  one  order  of  magnitude  smaller  than  the  first  term. 


In  view  of  other  uncertainties,  it  was  decided  to  drop  these  correcting  terms 
altogether,  simplifying  further  the  coi^putations.  A  similar  order-of- mag¬ 
nitude  analysis  indicated  that  both  <2  t  1-1) /4A^  and  10"®.  N  are  always  at 
least  one  and  three  orders  of  magnitude  smaller  than  unity  respectively, 
when  the  assumption  is  made  that  the  order  of  S  is  10“^.  Therefore,  we 
also  dropped  thca**  terms  from  thp  parenthesis,  with  the  final  result 

T  =  -10"®  .  N  .  A.  (F  (t  )-E.  F(t))  (38) 

o  o 

which  is  basically  similar  to  Rowlandson's  (ibid)  formulation  for  moderate 
to  long  ranges,  but  wMch  remains  accurate  at  short  ranges.  This  analysis 
also  shows  that  it  is  not  practical  to  improve  further  on  Freeman's  approach 
by  keeping  second-order  terms. 

A  Real-Time  Iterative  Solution  for  the  Bending 

Equation  (38)  still  suffers  from  the  disadvantage  of  Freeman's  original 
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ever,  if?  neither  necessary  nor  flesirabfe.  As  airrotlv  pointeci  out  by'’ {’ref- 
man  (il>id;  in  hia  original  paper,  using  the  radar  “ievation  aJigle  in  placf 
of  /3  tends  to  uriricrtostimate  T  because  the  coinpatiog  line  of  sight  (ends  to 
be  a*&ove  the  real  path.  Using  the  true  value  of  ^  ,  while  improving  the 
approximation,  tends  to  overestimate  the  bending  3.3  the  .slant  path  now  lies 
below  the  true  path.  We  .should  also  point  out  here  that  another  and  even 
larger  aoui'ce  of  error  comes  from  the  use  of  approximation  (26)  at  the 
longer  ranges.  The  real  elevation  error  e  =  and  the  true  bending  T 

are  related  by  the  expression 


_ _  j  COST  +  sir  T  .  tan  B  -  n/n^ 

arc  an  ;~lan“B  T’irn  t  -  cos  t  . 


tan  B 


(39) 


Since  T  is  always  a  small  angle,  if  we  make  the  cu-stomary  small-angle 
approximations  and  keep  up  to  second-order  terms,  using  (24),  (22),  (32) 
and  (37) 


T^/2 


T  .  tan  B  -  10 


N 


(1  -  E) 


<  = 


T  +  tan  B  -  tan  B 
o 


(40) 


Further,  Weisbrod  and  Anderson  (1959)  have  shown  that  over  short  ranges 

(41) 


-2  .  (N  -  N)  .  10“®  =  -lO'®  .  N  .  (1  -  E) 
o  o 


tan  B  +  tan  B 


tan  B 


2, 


so  for  short  ranges  the  numerator  of  (40)  reduces  to  r  72.  At  the  same 


time,  (tan  B  -tan  B)  «T,  so  the  denominator  of  (40)  tends  to  T. 
o 


t/2  (short  ranges) 


(42) 


and  this  approximation  will  undercorrect  longer  ranges,  since  as  the  range 
increases  indefinitely,  tan  B  increases  without  bounds  and  from  (40) 


T  (long  ranges) 


(43) 


If  we  apply  correction  (42)  to  the  initial  elevation  angle  B  only  once,  there 
will  be  a  tendency  for  compensation  at  long  ranges  since  f42)  undercompen¬ 
sates  at  long  ranges  which  is  exactly  what  we  need. 
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This  compensation  mechanism  indeed  works,  so  the  final  programming  was 
made  by  setting  -  B  initially  and  computing  an  approximate  bending 
by  using  (38).  The  initial  elevation  guess  is  corrected  by  adding  half  the 
approximate  bending  Tj^/2  and  the  final  bending  T  is  computed  using  (38) 
again. 

The  computer  requirements  are  very  small  and  the  only  subroutine  involved 
is  the  computation  of  the  single  exponential  E  in  (37).  Once  the  total  bending 
T  is  available,  all  the  other  parameters  can  be  extracted  from  it.  Accord¬ 
ingly,  Appendix  A  contains  the  calculations  for  Range  Error,  Elevation 
Error,  and  Doppler  Velocity  Error  Angle.  There  also  appear  suggestions 
for  quadratic  exponential  atmospheres,  and  for  a  possible  use  of  the  radar 
beam  as  an  atmospheric  probe.  We  regard  this  as  "spillout"  from  the  pre¬ 
sent  work. 

Comparisons  of  Retraction  Corrections  at  Low  Angles  of  Elevation  Using 
the  Proposed  Real-Time  Computations  with  Those  Obtained  by  Standard 
Solutions 

Tables  8  and  9  present  a  comparison  between  the  exact  solutions  for  tlie 
refraction  corrections  at  low  angles  of  elevation  for  a  sliUidard  CRPL  ref¬ 
erence  atmosphere  and  those  developed  in  this  report.  For  shortness,  we 
refer  to  the  exact  and  the  iteration  results. 

Tables  8  and  9  illustrate  the  remarkably  good  fit  between  the  iteration  and 
the  exact  solutions.  The  data  is  presented  in  tabular  form  beenuse  it  would 
not  have  been  practical  to  present  it  in  graphical  form,  as  the  two  lines  would 
have  plotted  one  on  top  of  the  other  in  most  cases.  The  fit  is  particularly 
good  at  medium  and  short  ranges,  an  area  where  Rowlandson's  (1968)  simp¬ 
lifications  start  losing  validity.  The  largest  absolute  error  in  range  error 
stays  under  3  meters  in  the  worst  case,  which  is  smaller  than  the  size  of 
most  targets.  Even  at  zero  elevation,  the  bending  and  elevation  error  dis¬ 
crepancies  do  not  exceed  .  2  mil,  which  are  smaller  than  differences  that 
occur  due  to  uncertainties  in  our  knowledge  of  the  state  of  the  atmosphere. 
This  is  especially  impressive  if  we  realize  that  the  iteration  technique  is 
perfectly  general  and  uses  no  "ad  hoc"  parameters  to  improve  the  fit. 

It  must  be  pointed  out.  however,  that  these  comparisons  were  made  with  a 
well-behaved  model,  and  that  greater  differences  may  be  expected  with  a 
real  atmosphere,  a  point  which  will  be  taken  up  later  when  the  discussion  of 
the  determination  of  the  decaj  constants  is  considered.  A  most  iimportant 
fact  is  that  only  about  107  multiplications  and  divisions,  70  additions  and 
subtractions,  2  exponential  subroutines,  2  square  roots  and  3  logical  de¬ 
cisions  arc  needed  for  a  complete  solution  of  all  refraction  errors,  which 
makes  the  procedure  near  real-time  even  on  ihe  smallest  computer.  The 
algorithm  is  much  faster  that  Rowlandson's  as  there  is  no  need  to  integrate 
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TABLE  8 


TABLE  9 


Comparison  between  elevation  and  Doppier  velocity  angle  error  for  a  standard 
CEPIatmosphere  for  different  ranges  and  radar  elevation  angles,  and  the 
values  obtained  by  the  iteration  real-time  technique 


Legend:  B  «  exact  value  of  error 

1  =3  error  computed  using  the  iteration  technique 
initial  radar  elevation  angle 
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numerically  for  the  elevation  error. 

The  Anisotropic  Atmosphere  -  A  Ray-Steering  Program 

The  detailed  refractometer  cro3S  sections  available  in  this  project  have 
shown  that  the  assumption  of  stratification  is  valid  only  as  a  first-order 
approximation.  Isopleths  of  refractive  index  (refractopleths)  are  generally 
tilted  with  respect  to  the  earth's  geopotential  surfaces  and  the  gradients 
are  far  from  isotropic  in  the  horizontal.  These  departures  exert  their 
greatest  influences  at  low  angles  of  incidence  for  the  propagated  rays  which 
are  of  interest  in  this  report. 

In  trying  to  follow  a  ray  through  an  anisotropic  medium,  we  are  confronted 
with  two  practical  problems.  The  first  one  is  the  description  of  the  medium 
itself  since  we  now  have  to  specify  the  three-dimensional  distribution  of  the 
refractive  index,  and  the  data  handling  problem  becomes  formidable.  The 
second  one  is  that  the  integration  of  the  ray  equations  has  to  be  accomplished 
by  some  sort  of  numerical  integration  scheme,  since  an  analytical  descrip¬ 
tion  of  the  complex  refractopleth  patterns  is  a  very  difficult  task.  Initially, 
we  considered  a  simplified  scheme  to  describe  the  observed  patterns  by  using 
orthogonal  polynomials,  following  the  lead  of  Wadsworth  (1951)  who  had  some 
uccoss  in  doing  the  same  with  pressure  patterns.  The  objective  in  mind 
was  to  relieve  the  storage  problem  in  the  computer.  In  common  with  what 
was  found  by  Lopez  and  Nason  (1967)  in  a  modelling  study  of  pi*ecipitation 
patterns,  which  resemble  refractopleths  a  great  deal,  this  approach  did  not 
prove  practical. 

A  description  of  the  patterns  by  gridpoints  is  the  next  approach,  and  would 
have  been  the  one  used  except  that  difficulties  in  procuring  an  adequate  com¬ 
puter  limited  severely  the  storage  capacity  available.  The  logical  compro¬ 
mise  was  to  use  the  human  brain,  wherein  a  hybrid  graphical -computational 
approach  was  used  in  the  exploration  work.  Two  independent  analyses  were 
performed  on  the  same  data  to  identify  human  error  in  pattern  display. 

All  distances  were  normalized  in  terms  of  the  radius  of  the  earth,  and  a 
practical  xinit,  the  "mer"  (for  milli-earth-radius),  equal  to  the  distance  sub¬ 
tended  by  one  mil  angle  at  the  earth's  surface  was  introduced  to  simplify  the 
computations.  A  simple  piecewise  integration  scheme  using  a  variable  hori¬ 
zontal  step  unit  was  next  devised  to  compute  the  terminal  elevation  of  the  ray 
at  the  end  of  each  step.  The  details  of  the  procedure  are  given  in  the  next 
paragraph. 

As  long  as  the  data  has  to  be  plotted  in  cross-section  form  to  check  its  valid¬ 
ity  in  exploration  work,  this  hybrid  procedure  is  perfectly  satisfactory.  It 
is  very  easy  to  computerize  the  whole  operation  by  describing  the  patterns 
by  means  of  grid  points,  storing  them  in  the  computer  and  letting  the  computer 
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do  the  differencing  and  interpolation. 

It  is  then  possible  to  envision  a  workable  system  in  which  two  airplanes  are 
stacked  one  sdsove  the  other,  one  porpoising  through  the  lower  atmosphere 
and  the  other  around  the  inversion,  both  telemetering  the  data  in  real-time 
to  a  computing  center  which  would  process  and  store  it  in  readiness  for 
tracing  any  selected  ray.  A  total  time  of  some  45  minutes  should  be  needed 
to  run  a  given  cross  section,  and  since  horizontal  differences  are  less  marked 
tfian  vertical  ones,  a  triangular  run  between  the  various  radar  centers  should 
define  the  horizontal  changes  in  a  practical  manner. 

Computations  for  Ray  Steering 

Computations  for  the  ray  tracings  were  carried  out  by  a  piece -wise  integra¬ 
tion  of  equation  (1)  as  follows.  An  arbitrary  interval  AHD^  in  mers  was 
selected  at  the  ith  point,  so  that  if  the  ray  were  extrapolated  by  this  much, 
the  vertical  gradient  remained  sensibly  constant.  It  was  then  assumed  that 
1)  the  difference  between  the  ith  ray  segment  and  its  horizontal  projection 
was  a  second  order  quantity.  This  assumption  is  quite  tenable  for  rays  up 
to  at  least  100  mils  elevation;  2)  the  ray  can  be  approximated  by  an  arc  of 
circle  whose  curvature  equals  the  average  curvature  of  the  segment. 

The  curvature  at  any  point  in  the  segment  equals  the  local  logarithmic  deriv¬ 
ative  of  the  refractive  index  (l/n)  .  (dn/dR  ),  where  R  is  the  local  radius 
of  curvature.  If  we  substitute  the  averaged  derivative  over  the  segment  by  the 
corresponding  finite  difference  equivalent,  noting  that  the  factor  (l/n)  differs 
from  unity  by  at  mos^  third-order  quantity  and  that  n  =  1  +  10"?  N,  we  ob¬ 
tain  the  expression  An.  10“6/AR  for  the  average  curvature.  Refractive 
indices  are  available  only  in  the  form  of  cross-sections  with  limited  data, 
so  a  compromise  is  needed  for  the  selection  of  the  radius  of  curvature  incre¬ 
ment  AR^.  If  we  make  it  small  to  approximate  the  local  derivatives  we  will 
have  uncertainties  in  the  determination  of  the  averaged  refractive  index  change 
AN.  If  we  make  it  too  large,  we  will  smooth  out  many  real  irregularities. 

In  the  end,  a  value  of  AR  =  .  01  mers  (about  67  meters)  was  chosen  as  the 
best  compromise.  For  tfiis  fixed  value  of  AR  ,  the  average  curvature  for 
t^  segment  AHD  is  10"^  .  AN  and  the  average  local  bending  will  be  At  = 

AN  .  AHD/ 10  (mils). 

The  net  change  in  elevation  from  the  extrapolated  position  then  follows 
immediately  since  it  equals  the  difference  in  elevation  due  to  the  opposing 
effects  of  the  earth's  and  the  ray's  curvature.  Since  both  trajectories  are 
assumed  circular,  the  computations  are  greatly  simplified.  The  angle  of 
emergence  follows  from  relatio.n  (11)  between  the  bending  and  accumulated 
horizontal  distance. 

The  scheme  of  the  computations  for  the  ith  point  is  as  follows:  (All  distances 
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Computations  are  iterated  until  the  nth  and  last  point,  and  then  the  refraction 
corrections  are  obtained  from 

RGERR  =  10"^ 

i=0 

T  =  HD  +  B  -  B  (mils)  (45) 

non 

€  =  B  +  HD  /2  -  10^  .  h  /HD  (mils) 

on  n  n 

The  range  error  "RGEER"  can  be  changed  to  more  conventional  units  by  using 
the  equivalences  1  mer  6370  m  20900  ft.  As  a  check,  the  value  obtained 
for  the  total  bending  from  (45)  must  coincide  with  the  value  of  in  (44). 

The  computations  involved  in  (44)  are  extremely  simple,  and  by  exercising 
judgement  over  the  selected  intervals,  we  found  that  rays  could  be  traced  very 
rapidly  with  as  few  as  a  dozen  key  points  along  the  trajectory.  We  have  checked 
the  approximations  involved  in  (44)  by  repeating  the  calculations  using  shorter 
intervals  and  more  accurate  formula,  and  found  that  as  long  as  the  initial  ele¬ 
vation  is  under  100  mils  there  is  no  need  for  greater  sophistication.  Currently 
available  inputs  just  do  not  justify  the  use  of  greater  computing  effort. 


SECTION  VI 


SUMMAKif  AND  RECOrvIMEI^DATIOMS 


Summary 

With  prevailing  northwest  wind,  the  inversion  immediately  to  the  south  of 
Tranquillon  Peak  is  likely  to  be  destroyed  by  the  hydraulic- jump  effect 
associated  with  sufficiently  high  Froude  numbers.  This  will  often  be  noticed 
as  an  erosion  of  the  cloud  deck  at  inversion  height.  This  circumstance 
should  allow  excellent  target  tracking.  The  well-mixed  atmosphere  here 
will  approximate  a  standard  atmosphere.  The  mixed  condition  should  exist 
some  20  miles  to  seaward. 

Generally,  the  area  immediately  to  the  west  of  Tranquillon  Peak  will  have 
low-level  mixing  due  to  eddying  around  the  land  mass,  in  which  case  one  will 
find  a  clear  sky  area  along  the  beach  where  the  mixing  has  eroded  the  cloud 
deck  (present  at  inversion  height).  At  higher  altitude  one  will  find  a  set  of 
standing  waves  from  the  disturbance  of  the  pesds  itself. 

In  the  northerly  sections  wave  trains  across  Santa  Inez  Valley  can  be  expected. 
The  inversion  may  also  have  a  warped  geometry  which  seemingly  climbs  up¬ 
ward  to  meet  Tranquillon  Pesdt,  This  warpir.  {  can  actually  be  seen  by  a 
ground  observer  on  the  north  hills  when  the  ’  ow  morning  sun  scatters  off  the  . 
haze  layer.  The  trapping  and  bending  significance  of  a  warped  or  tilted  inver¬ 
sion  layer  has  been  brought  out  in  the  section  on  radio  hole  missions. 

There  will  be  a  slight  elevation  difference  in  the  inversion  between  the  north 
and  south  sides  of  San  Nicholas  Island.  The  gradient  laj’er  will  bulge  over  the 
island  on  the  order  of  1,  000  ft.  It  will  have  a  V-wake  structure  propagating 
out  to  the  north  and  soudi  sides  of  the  Island.  It  will  have  a  large  standing 
wave  off  the  windward  nose  of  the  island.  There  is  strong  convergence  of  the 
gradient  layer  eastward  of  the  mid  point  of  the  island  and  a  hydraulic  jump 
induced  mixing  of  the  Inversion  layer  downwind  of  this,  producing  a  well  mixed 
wake  off  the  eastern  end  of  the  Island. 

The  inversion  at  Point  Mugu  may  be  washed  out  from  eddy  mixing  associated 
with  shear  where  the  sea  breeze  front  meets  the  north  westerlies.  It  can  also 
be  multi  layered  or  lensed  under  moderately  weak  stability  conditions  (to  be 
seen  in  Figure  16,  Section  IV,  8  June  example).  The  seaward  strength  of  the 
inversion  may  be  experienced  close  along  shore  and  over  the  station  at  Point 
Mugu  when  strong  stability  conditions  are  present. 

The  statistics  of  refractive  index  corrections  so  far  accumulated  suggest  that 
corrections  increase  from  shore  to  seaward  as  a  function  of  increasing  sharpness 
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models  are  employed  produce  greater  corrections  tlian  those  produced  from 
Raob  or  aircraft  spiral  data  employing  "stratified  atmosphere"  models.  The 
average  separation  of  relative  error  between  the  extremes  of  data  and  form¬ 
ulae  employed  in  a  comparison  cf  refractive  corrections  between  Raobs. 
spirals,  and  cross  sections  is  68  percent  over  six  cases  analyzed  in  detail. 
The  greatest  differences  are  found  in  the  elevation  errors. 


Investigation  of  the  relation  between  the  altitude  position  of  a  C>band  beacon 
equipped  aircraft  and  its  indicated  radar  position  when  loss  of  lock  is  exper¬ 
ienced  shows  that  Uie  actual  position  of  signal  loss  is  much  better  defined  if 
a  full  cross  section  and  steering  techniques  are  used  to  correct  the  radar  for 
refraction. 

Signal  lose  is  associated  with  atmospheric  wave  nodes  and  sna^'ply  converging 
refractopleths  in  the  gradient  layer.  This  accounts  for  crossing  over  of  an¬ 
tenna  angles  in  the  vertical  plane.  This  is  reflected  in  the  radar's  hunting  in 
the  vertical  plane  at  very  low  angles  (<  10  mils).  This  does  not  preclude  that 
hunting  may  also  be  a  function  of  the  usually  alleged  multipath  (which  in  these 
circumstances  is  almost  synonomous). 

The  atmospheric  waves  found  at  sea  at  the  inversion  interface  may  dwell  for 
considerable  time  in  one  place  (8C  min. ).  Where  associated  with  the  Tranq- 
uillon  Peak  topographically  induced  waves  the  node  positions  of  the  waves  may 
shift  in  the  lateral  plane,  but  wave  size  and  geometry  remains  the  same  (45 
min.  comparisons). 


Examination  of  REFRAC  correction  program  supplied  by  Western  Test  Range 
reveals  that  the  program  is  mathematically  correct  but  overly  complex.  Its 
inability  to  produce  a  satisfactory  run  forced  a  redesign  of  the  program  which 
now  contains  only  41  computing  statements:  involving  77  multiplications  and 
divisions,  60  additions  and  subtractions,  ?  square  root  routines,  and  11  logical 
decisions  per  layer.  This  is  the  sole  requirement  for  the  real  time  correction 
of  a  non-trapped  ray.  The  program  has  run  on  a  TELCO MP  system  serviced 
by  a  POP- 7  computer.  The  program  ser  ves  any  elevation  angle,  including 
negative  angles  and  warns  of  impending  trapping  conditions. 

The  problem  of  steering  rays  through  an  anisotropic  atmosphere  is  described 
and  formulae  presented  for  appropriate  ray  tracing  where  strong  irregularities 
exist  along  the  sighting  line. 

Recommendations 


If  accuracy  and  reliability  demands  on  the  trilateration  system  at  Western  Test 
Range  become  no  greater  than  they  presently  are,  do  nothing  except  learn  more 
about  the  hazards  of  side  lobe  splitting  as  a  function  of  the  inversion  layer. 
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This  can,  and  may  have  caused  at  least  one  missile  destruction  in  the  past. 

Use  a  mid-path  reiiector  van  auernallve  steeper  angle  path)  for  checking 
beacon  fade  on  the  pre- launch  Minuteman  checks  when  inversion  interfer¬ 
ence  iu  suspected* 

If  accuracy  and  reliability  demands  on  the  trilateration  system  at  Western 
Test  Range  are  likely  to  become  greater  in  the  near  future,  then 

1)  Confirm  the  inversion  effects  brought  to  light  In  this  investigation. 

2)  Introduce  meteorological  forecasts  of  the  condition  of  the  inversion 
before  shot  time  (see  Appendix  B). 

3)  Forget  about  the  tiny  differences  in  the  various  correction  formulae 
and  devise  better  techniques  for  coping  with  the  very  real  differences  in  at¬ 
mospheric  data  inputs. 

4)  Go  to  real-time  refraction  corrections  with  built  -in  trapping  warn¬ 
ings. 

5)  Focus  attention  on  procuring  full  cross  sections,  starting  30-40 
minutes  before  shot  time. 

We  propose  for  pre-launch  cross  sections  that 

1)  Two  aircraft  begin  exploration  of  the  inversion  at  sea,  at  t-S^O  (t  = 

0  =  shot  time)  where  low-angle  rays  will  intercept  the  gradient  layer  (depends 
on  the  height  of  the  layer).  The  aircraft  will  establish  layer  height,  gradient, 
and  geometrical  irregularity. 

2)  At  t-35  one  aircraft  begins  to  map  a  cross  section  for  the  first  15 
miles  in  front  of  the  radar,  from  the  top  of  the  inversion  to  the  surface.  The 
second  maps  the  precomputed  interception  zone  (where  the  rays  will  pass 
through  the  inversion). 

3)  During  the  entire  pre -shot  mission,  both  aircraft  telemeter  their 
data  directly  to  the  beach  in  real  time.  The  cross  sections  are  constructed 
in  computer  memory.  Prior  to  the  t-35  mappings,  preliminary  ray  tracings 
are  made  to  help  identify  regions  for  concentration  of  mapping  effort  during 
the  next  30  minutes.  At  t-5  the  mapping  aircraft  are  diverted  off  the  intended 
trajectory.  Between  t-5  and  t  =  0  the  cross  sections  are  ray-traced  again 
for  possible  danger  points.  At  t  =  0  the  cross  sections  are  used  for  real¬ 
time  refraction  correcticwis  to  missile  trajectory. 

Appropriate  Exploratory  Work 

It  has  been  learned  that  radar  targets  have  been  placed  on  Santa  Cruz  Island 
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for  calibration  in  the  past  (Bankston,  1964).  Such  could  also  be  useful  at 
San  Nicolas  Island  or  at  sea  with  a  balloon  and  spar  buoy  combination, 
where  greater  height  is  required,  tethered  balloons  monitored  for  position 
by  cine  theodolites  should  be  considered.  These  fixed  targets  would  then  be 
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sections,  management  of  aircraft,  and  instrumentation  and  telemetry,  and 
ray  tracing  for  error  prediction.  We  further  propose  a  system  be  investi¬ 
gated  wherein  a  backward  solution  of  the  refraction  equations  is  performed 
so  thv..  .Adjusting  parametei  s  are  computed  from  the  cross-section  data. 

This  approach  has  developed  as  unexpected  spill-out  from  the  present  project. 
We  have  devised  a  preliminary  method  for  performing  the  mathematics  of 
this  back  solution  which  requires  only  knowledge  of  the  true  position  of  two 
targets  at  two  different  ranges.  'The  radars  themselves  could  then  act  as 
atmospheric  probes.  During  the  research  aspects  of  such  a  program,  two 
or  three  high-speed  horizon  level  and  near-horizon  fly-bys  should  be  incor¬ 
porated  wherein  the  tracking  of  a  true  inertial  target  reference  should  reveal 
anomalous  sections  of  the  atmosphere,  manifested  as  unstable  radar,  to 
which  would  be  applied  backward  solutions. 


Summary  Statement 


The  era  of  mathematical  exercisefj  in  designing  correction  models  based  on 
simple  assumptions  is  over.  When  two  or  three  mil  errors  can  be  demon¬ 
strated  solely  due  to  atmospheric  anomalies,  it  makes  little  difference  that 
mathematical  model  A  is  better  than  mathematical  model  B  because  it  com¬ 
putes  to  .  01  mils.  If  greater  accuracies  are  to  be  expected  at  low  angles 
from  the  Western  Test  Range,  we  must  recognize  the  fallacy  of  assuming 
an  ideal  standard  atmosphere. 
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All  r«fraotlon  errors  At**  interpendent.  For  the  iteration  technique,  we  have 
selected  the  bending  error  as  the  base  from  which  all  tiie  other  refraction 
errors  are  obtained,  since  the  convergence  of  the  iterations  for  the  bending 
error  is  very  fast  and  the  computing  procedure  straightforward.  In  what 
follows,  we  show  in  detail  the  computation  of  the  remaining  refraction  errors. 

RANGE  ERROR 

Initially,  we  neglect  the  curvature  of  the  path  in  the  range  computations 
(Section  V).  The  elementary  retardation  error  is  (Weisbrod  and  Anderson, 
1959) 

d(RGERR)  =  -(n  -l)dS  (46) 


Introducing  (24),  (22),  (32)  and  (33) 


d(RGERR)  -  - 


N  .  exp  (t^  -  t^)  .  dt 
o  o _ 

~AT~C03~BZ 


which  can  be  integrated  directly  to  yield 

RGERR  =  -10‘®  .  N  .  sec  3  .  (F(t„)  -  E  .  F  (t))/2A^ 
o  o  o 

where  F(x)  is  defined  by  (36). 

Comparing  (48)  and  (38)  we  obtain 

RGERR  -  T  .  sec 

If  we  now  add  the  correction  for  curvature  of  path  (17) 

RGERR  =  (T  .  sec  +  S  .  T^/24)  .  r^ 


(49a) 


where  the  range  error  has  been  multiplied  by  to  expre.ss  it  in  conventional 
units.  We  can  see  from  (49a)  that  for  small  angles  of  elevation  the  choice  of 
3  is  not  very  critical  since  the  secant  varies  very  slowly  at  low  angles  of 
ePevation.  We  can  therefore  substitute  for  3^  in  (49a)  and  obtain  as  the 
final  expression 

RGERR  =  T  .  r  .  (sec  B  /2A^  +  S  .  t/24)  (50) 

o  o 
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ELEVATION  ERROR 


An  iterative  procedure  can  be  useW  to  determine  the  elevation  correction  "t". 
From  equation  (40)  we  can  aee  that  once  T  is  available,  the  elevation  correc¬ 
tion  c  depends  only  on  an  estimation  of  the  emergence  angle  B.  Although 
(25)  and  (28)  provide  the  basis  for  an  initial  estimation  by  setting 
the  angle  B  is  generally  small  and  is  not  well  determined  by  its  cosine.  It 
is  better  to  modify  the  computations  in  a  manner  similar  to  that  done  in 
Section  V  and  to  proceed  as  follows: 


Let 


G  =  (tj  -  t^)/2A^  -  10"®  .  N  .  (1  -  E) 


where  t,  t  and  E  are  obtained  from  (31)  and  (37)  by  setting  ^  =  B 
O  ^  o  o 


and  let 


K1  =  1  -  cos  B 


(51) 


(52) 


be  computed  accurately  by  using  a  suitable  rational  approximation  as  suggested 
in  Section  V. 


Then 

KB  =  BI  +  G  .  (1  -  KI) 

(53) 

and 

tan  B  =  (KB  . 

(2  -  KB))^^^/(1 

-KB) 

(54) 

tan  B  =  (KI . 

(2  -  KI))^^^/(1 

-  KI) 

(55) 

A  first  estimation  is  now  obtained  by  using  (40).  A  second  iteration  is 

(56) 


now  made  by  recomputing  t,  t^  and  E  from  (31)  and  (37)  using 


/3  =  B  +  €, 
o  o  1 


This  second  iteration  yields  a  final  estimation  for  the  elevation  error  c. 

The  convergence  of  this  Iterative  procedure  is  sufficiently  fast  that  no  further 
meaningful  improvement  in  €  can  be  expected  from  further  iterations. 

DOPPLER  VELOCITY  ERROR  ANGLE 


Miilman  (1961)  has  shown  that  doppler  velocity  error  in  the  radial  direction 
which  can  be  attributed  to  refraction  is,  to  a  first  order  approximation, 
measured  by  the  angle  (5  (Figure  21),  which  is  the  angle  between  the  ray 
direction  and  the  slant  path.  From  Figure  21,  it  can  be  seen  that 


6  =6 


T  -  B 


(57) 
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(58) 


f 


r 


I--- 


But  since  -  B  =  € .  the  elevation  error,  we  have  that 
o  o 

6  =  C  -  T 


which  shows  that  a  doppler  correction  angle  is  immediately  available  from 
the  estimates  for  the  total  bending  and  me  elevation  tmgle  error. 

A  QUADRATIC  EXPONENTIAL  MODEL 


In  our  previous  work,  an  algorithm  has  been  developed  to  perform  fast  and 
accurate  refraction  corrections  for  an  atmosphere  whose  refracting  index 
decays  exponentially  with  heij^t.  The  logical  question  has  been  considered 
by  Bean  and  Thayer  (1959)  and  practical  experience  has  shown  that  the  ex¬ 
ponential  model  fits  most  observed  profiles  reasonably  well.  Rowlandson 
U968)  has  shown  that  even  large  departures  from  the  exponentisd  model  re¬ 
sults  in  modest  differences  in  the  computed  corrections  if  suitable  parameters 
are  chosen. 

We  can  use  available  direct  inputs  such  as  radiosonde  data,  spirals,  etc.  , 
to  derive  the  parameters  of  a  best-fitting  exponential  model.  There  is  an¬ 
other  alternative  (suggested  in  our  summary  and  recommendation)  if  targets 
of  known  position  are  available,  of  using  the  radar  as  an  atmospheric  probe 
by  performing  a  back  solution  of  equation  (38).  For  the  moment,  however, 
we  are  concerned  with  the  first  alternative. 

If  a  refractive  profile  is  available,  the  solution  for  the  best  parameters  is 
handled  by  standard  curve-fitting  techniques  such  as  a  least-squares  routine. 
The  simple  exponential  model  which  we  have  been  considering  involves 
fitting  a  least-squares  straight  line  to'H^  and  ln(N. ).  As  soon  as  several  of 
these  plots  are  made,  however,  one  is  immediately  tempted  to  improve  the 
fit  by  going  to  higher-order  polynomials,  especially  where  a  strong  inversion 
is  present.  If  the  monomial  in  the  exponent  of  equation  (22)  is  substituted 
by  a  binomial,  we  obtain  a  quadratic  model  where 

N  =  .  exp  (-Dj  .  H  -  D2  .  H^)  (59) 

and  we  could  fit  a  parabola  to  the  plot  of  H^,  In(N^)  instead  of  a  strai^t  line. 

If  we  make  the  appropriate  substitutions  in  (22)  and  (25)  and  neglect  terms 
higher  than  second-order  in  the  subsequent  derivations,  we  find  that  we  only 
need  to  make 

A  =  (Dj/2)^^^  .  (1  +  Dg  .  tan^  0^/(Dj/2))"^^^  (60) 

instead  of  the  value  given  in  (31)  to  arrive  at  the  same  solution  (38)  for  T. 

This  makes  the  use  of  a  quadratic  model  almost  as  practical  as  that  of  the 
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simpler  linear  exponential,  since  the  fitting  of  a  parabola  to  a  set  of  data  is 
still  a  '/ery  simple  nomputational  procedui^.  However,  in  general,  we  have 
seen  that  single  point  radiosonde  data  may  not  detect  the  atmospheric  anum- 
aUes  which  so  much  affect  tracking  at  low  angles  of  elevation,  so  the  gain 
in  accuracy  using  a  quadratic  model  may  be  illusory.  If  detailed  cross- 
section  data  or  targets  of  known  position  are  available,  it  seems  preferable 
to  use  die  linear  exponential  model  to  estimate  an  equivalent  average  decay 
pareuneter  in  the  manner  suggested  in  the  next  paragraph. 

DETERMINATION  OP  THE  DECAY  CONSTANTS  - 
POSSIBLE  USE  OF  THE  RADAR  BEAM  AS  AN  ATMOSPHERIC  PROBE 

It  is  obvious  that  the  success  of  the  real-time  computations  depend  on  a 
proper  estimation  of  the  parameter  N  and  the  decay  constant  DEN  in  the 
case  of  a  simple  exponential  atmosphere  or  the  two  decay  constants  D^  smd 
D-  of  the  quadratic  model.  A  least- squares  fitting  routine  could  be  used  to 
d^ermine  the  constants  when  Inputs  such  as  radiosonde  data  are  available. 
Rowlandson  <1968)  has  discussed  this  problem  in  more  detail. 

In  the  course  of  our  work,  we  had  atmospheric  data  at  our  disposal,  and  it 
soon  became  evident  that  the  radar  inputs  themselves  could  be  used  to  esti¬ 
mate  the  decay  constants  if  we  could  solve  the  refraction  equations  back¬ 
wards.  We  posed  the  question:  If  we  know  what  the  elevation  angle  and 
bending  corrections  should  be,  what  average  state  of  a  hypothetical  stratified 
atmosphere  would  have  such  a  solution  for  the  given  initial  conditions  ?  Since 
the  equations  are  non-linear  and  the  unknowns  are  implicit,  some  sort  of 
Iterative  solution  is  necessary.  We  have  made  a  preliminary  exploration  of 
the  problem  and  have  succeeded  in  devising  a  technique  based  on  the  bisection 
method  to  Isolate  the  roots  rather  rapidly.  The  procedure  appears  to  be 
both  convergent  and  stable. 

The  bonus  to  be  gained  from  such  a  procedure  is  quite  apparent.  If  two  tar¬ 
gets  of  known  position  are  available,  radar  inputs  define  the  constants  of  an 
average  equivalent  atmosphere.  The  targets  can  be  fixed  or  stationary,  but 
in  view  of  the  many  desirable  look  angles,  a  fly-by  using  inertial  guidance 
seems  to  offer  the  greatest  potential  to  define  an  equivalent  atmosphere. 

We  have  already  made  computations  using  the  known  position  of  the  aircraft 
as  defined  by  the  radio  altimeter  and  the  radar  range  as  the  best  available 
absolutes,  and  have  come  out  with  decay  constants  which  can  reproduce  the 
actual  departures  at  intermediate  points  far  better  than  the  radiosonde  of 
simple  spiral  data.  Since  the  procedure  cannot  pinpoint  the  exact  location 
of  the  anomalies,  but  averages  them  throughout,  this  short  survey  indicates 
that  the  observations  should  be  carried  at  two  elevations  and  two  ranges,  one 
at  least  a  fairly  long  one  to  insure  that  a  thick  slice  of  the  lower  atmosphere 
is  traversed.  This  is  an  area  that  should  be  explored  further. 
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A.PPENDIX  B 


Pacific  inversion  layer  activity 


During  our  1967  Vandenberg  operaticms  serious  thought  was  given  to  the 
possibility  of  predicting  inversion  layer  activity  (changes  in  strength  and/or 
geometry)  by  statistical  techniques  and  mathematical  modeling.  Chie  would 
hope  to  identify  parameters  for  which  physical  observations  and  theory  would 
in^cate  they  had  strong  control  over  inversion  activity.  The  control  wotild 
be  confirmed  by  statistical  correlation  from  which  one  would  next  model  the 
theoretical  association  of  these  parameters.  The  intent  would  be  to  simulate 
inversion  activity.  The  applied  reason  for  such  effort  would  be  for  predicting 
inversion  conditions  that  were  hazardous  to  missile  tracking.  The  first 
phase  of  this  progrsnn  has  been  completed,  that  of  identifying  parameters  for 
which  theory  or  past  work  indicated  a  degree  of  control  over  inversion 
activity.  A.  condensation  of  this  research  appears  in  this  Appendix.  A  good 
understanding  of  inversion  activity  and  its  significance  to  radar  can  be  obtained 
from  this  review. 

GENERAL  INVERSION  DESCRIPTION 

Tracking  variations  at  Vandenberg  Air  Force  Base  can  be  traced  to  the 
presence  of  the  strong  northeast  Pacific  inversion  through  which  most  tracking 
must  take  place.  The  inversion  layer  is  formed  at  the  interface  between 
subsiding  ?nti- cyclonic  air  from  aloft  and  the  marine  layer  generated  immedi¬ 
ately  above  the  ocean  surface.  These  layers  have  markedly  different  temp¬ 
erature  and  water  vapor  characteristics.  Sandwiched  between  these  two  air 
masses,  is  the  inversion  layer,  about  150  m.  in  thickness,  and  rising  from 
an  altitude  of  some  400  m.  off  the  coast  of  California  to  2,  000  m.  over 
Hawaii.  It  also  experiences  a  slight  upward  warping  along  the  immediate 
California  coast. 

Figure  B-1,  after  Neiburger,  1965,  shows  the  general  summer  position  of  an 
amt i- cyclonic  high  pressure  system  which  is  characterized  by  dry  subsiding 
aiir.  Note  the  streamlines  of  flow,  diverging  from  the  center  of  the  high 
pressure  area,  flowing  along  the  coast  of  California,  then  to  the  southwest  to 
form  the  northeast  trades. 

A  large  lateral  plane  of  demarcation  (the  inversion  layer)  between  the  upper 
subsiding  air  and  lower  marine  air  mass  is  maintained  by  the  temperature 
inversion  at  their  interface.  Since  the  temperature  increases  markedly  with 
height  in  this  layer,  strong  dynamic  stability  is  created.  This  stability  is 
generally  sufficient  to  prevent  penetration  by  the  convective  cells  generated 
by  the  heat  released  from  the  ocean  surface.  As  one  moves  westward  from 
the  coast  of  California,  however,  cumulus  convection  increases  due  to 
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gur*  B-l  NORMAL  THREE-DIMENSIONAL  TRAJECTORIES  OF  AIR  PARCELS  STARTING  AT  VARIOUS 
POSITIONS  AT  TOOmb.  (NUMBER  ALONG  TRAJECTORIES;  HEIGHT  IN  THOUSANDS  OF 
FEET  ).  AFTER  NEIBUR6ER,  1965 


. . . . . 


increasing  ocean  lemperaiure.  It  increases  to  a  degree  where  the  upward  trans¬ 
port  of  vapor  Euid  momentum  are  sufficent  to  penetrate  the  stable  layer.  The 
nffn«tr*>.tion  starts  with  mild  perturbations  of  the  layer  and  increases  to 
actual  penetratiem  of  the  layer  in  the  vicinity  of  Hawaii.  We  recog^ze  these 
penetratioris  as  chimney-type  ctunulus  clouds.  It  is  the  upward  transport  of 
vapor  and  momentum,  through  convection,  that  is  responsible  for  both  the 
raising  of  the  inversion  layer  as  one  goes  westward  and  the  eventual  destruc¬ 
tion  of  the  inversion  layer.  Since  greater  and  greater  aunounts  of  vapor  are 
mixed  into  the  dry  air  above,  the  dry  air  becomes  more  like  the  air  of  the 
lower  marine  layer.  The  inverslcm  layer  itself  gets  thicker  with  more  diffuse 
botmdaries  as  the  mixing  reaches  upward.  Eventually  the  inversion  demarcation 
completely  disappears.  The  interface  between  the  dry  and  moist  air  is  not 
a  discontinuity,  for  mass  transfer  actually  takes  place  between  the  dry  and 
moist  air.  It  is  the  physical  characteristics  of  the  mass  that  change. 

The  preceding  paragraphs  (taken  largely  from  Riehl.  1951.  and  Neiburger. 

1965)  should  give  one  a  gross  concept  of  the  nature  of  the  northeast  Pacific 
inversion.  This  is  the  crux  of  our  concern  with  respect  to  tracking  problems. 

The  balance  of  this  exposition  reviews  the  real  effects  of  the  inversion  as 
demonstrated  by  field  experiements.  by  Carr  et  al  (1958)  and  Bowen  et  al  (1965). 
and  supplements  this  with  statistical  compilations  by  these  same  authors,  and 
also  Sheddy  (1966).  Next  the  work  of  Riehl  et  al  (1951),  and  Neiburger  (1965) 
is  reviewed  in  sufficient  detail  to  isolate  the  physical  parameters  which  their 
data  indicate  could  reasonably  be  suspected  of  exercising  some  control  over 
the  activity  of  the  inversion  layer.  These  parameters  are  to  be  investigated 
for  their  value  as  predictors  of  inversion  layer  activity.  These  parameters 
are  tabulated  in  the  summary  at  the  end  of  this  writing.  An  investigative 
approach  to  the  handling  of  these  parameters  is  suggested,  along  the  lines  of 
an  eigenvector  analysis  by  Stidd  (1967);  or,  a  screened  multiple -discriminate 
analysis  by  Enger  (1965).  The  following  review  will  establish  the  state  of  the 
art  regarding  work  on  the  Pacific  inversion  pertinent  to  existing  tracking 
problems  at  Vandenberg  as  well  as  establish  a  point  of  departure  for  an 
investigation  regarding  inversion  activity  prediction. 

STATISTICAL  COMPILATIONS  AND  FIELD  EXPERIMENTS 
ON  REAL  EFFECTS  OF  THE  INVERSION 

As  an  indication  of  the  occurrence  of  strong  refraction  conditions  over  San 
Nicolas  Island  and  Point  Muga,  Carr  has  plotted  data  from  1956,  Figure  B-2 
showing  the  time  variation  (annual)  of  the  frequency  of  occurrence  of  inversions 
strong  enough  to  produce  trapping  conditions  (v/ith  refractive  index  greater 
than  36  N  per  1, 000  feet).  The  circled  numbers  adjacent  to  data  points  are 
numbers  of  observations.  Note  that  the  frequency  of  occurrence  of  tracking 
condtions  has  a  maximum  in  summer.  These  data  are  in  agreement  with 
similar  compilations  by  Sheddy  (1965)  and  Bowen  (1966),  among  others. 
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Additional  data  compiled  on  the  distribution  of  elevated  layers  during  the 
summer  maxima  are  shown  in  Figure  B-*3.  The  modes  of  these  distributions 
suggest  that  we  are  likely  to  find  a  layer  with  a  40  or  50  N  unit  change  with 
250  to  500  feet  of  layer  thiCfOiesa»  and  with  the  base  or  uie  rayex*  ii  uuu  icet 
In  altitude. 


In  the  work  of  Bowen  et  al  (1965).  we  find  £m  empirical  display  of  the  relation 
between  the  inversion  layer  and  trapping  conditions.  Bowen  and  his  colleagues 
set  up  three  receivers  on  the  abruptly  precipitation  coast  of  Point  Buchon 
170  miles  northwest  of  Los  Angeles.  These  stations  monitored  signal  strength 
from  Hawaiism  television  stations  at  several  frequencies.  In  Figure  B-4,  we 
note  from  the  signal  strength  record  (left  margin)  that  stations  A  and  C 
received  stronger  signals  the  closer  the  inversion  layer  came  to  their  res¬ 
pective  altitude  levels.  As  the  inversion  layer  climbed  with  time  the  strength 
of  the  signal  dropped.  These  experimental  recordings  accentuate  the 
importance  to  be  attached  to  inversion  layer  and  receiving  station  geometry. 

The  fluctuations  in  the  signal  strength  curves  are  alleged  to  result  from  the 
passage  of  undulations  in  the  geometry  of  the  inversion  layer  itself.  In 
general,  the  signals  that  were  received  at  the  upper  and  lower  sites  were  similar 
in.  fading  characteristics  but  not  in  power  levels  received.  Signal  level  at 
the  upper  site  was  higher  than  the  lower  site.  On  several  occasions  signals 
were  received  at  respective  stations,  only  when  the  inversion  layer  moved 
directly  through  the  antenna  structure  (ducting).  On  these  occasions  the 
inversion  layer  was  characterized  by  very  sharp  gradients. 


From  the  same  work  the  investigators  have  compiled  the  statistics  of 
signal  strength  received  as  a  function  of  time,  diurnally  and  seasonally.  We 
note  in  Figure  B-5  the  signal  strength  peaking  at  about  1530  in  the  afternoon. 

This  is  probably  near  the  time  of  maximum  absorption  of  insolation 
(temperature)  at  the  top  of  the  marine  layer,  and  conceivably  could  be  a 
contributing  mechanism  to  diurnal  fluctuation  of  the  layer  (this  parameter 
should  be  included  in  any  inversion  predictor  investigation), 

Sheddy  (1966)  has  made  a  study  of  the  magnitude  of  the  angle  of  arrival 
eiTor  of  a  horizontal  ray  using  data  from  radiosondes  from  different 
geographical  locations  (San  Diego,  Weather  Ship  N  in  mid-ocean  and  Hilo, 

Hawaii)  and  time  of  year.  The  angle  errors  computed  from  radiosonde  pro¬ 
files  are  then  compared  against  angle  errors  calculated  from  the  use  of  the 
surface  index  of  refraction,  in  the  CRPL  reference  atmosphere  model. 

Sheddy  shows  for  San  Diego  that  angular  error  is  greatest  in  summer  as 
might  be  expected  from  the  higher  water  contents  in  the  lower  atmosphere  during 
this  season;  and  Figure  B-6,  after  Sheddy,  shows  the  difference  in  error, 
e^,  between  a  radiosonde -based  calculation  and  CRPL  reference  atmosphere  - 
based  calculation  also  is  greatest  in  the  summer.  (Magnitude  of  error  is  lower 
at  higher  elevation  angles,  but  the  shape  of  curve  is  the  same.  There  is  also 
a  v/inter  night  anomaly  at  San  Diego,  where  the  error  difference  becomes 
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Fiflur#  B-3  DISTRIBUTION  OF  ELEVATED  LAYERS  WITH  THE  REFI 
GRADIENT  GREATER  THAN  ! -36N/I,000  FEET  I  DURING 
OF  I95S-57  OVER  SAN  Ni COLAS  ISLAND  (ll0  08SeR\ 


Figur*  B-5.  AVERAGE  DIURNAL  DUTY  CYCLE  FOR  CHANNEL  9.  AFTER  BOWEN,  1965. 


large,  possibly  due  to  strong  radiational  cooling  at  the  surface  producing 
a  meaningless  surface  N).  Note,  however,  that  the  more  seasonally  homo¬ 
geneous  atmosphere  at  Hilo,  Figure  B-7  has  a  correspondingly  flat  error 
distribuiiun;  and  has  a  generally  lower  magniluue  of  errux*  due  lu  iiie 
weaker  inversion-associated  water  vapor  gradient.  In  the  discussion  of 
Kiehl's  work,  peneti’ative  convection  (increasing  westward  towards  Hawadi) 
is  the  mechanism  for  transporting  moist  marine  air  through  the  inversion 
interface  to  the  dry  air  aloft.  This  process  tends  to  destroy  sharp  gradients. 
From  Sheddy's  graphs,  one  concludes  that  any  atmosphere  or  any  season  has 
acceptable  error  at  HUo  where  weak  gradients  exist.  The  mid -ocean  weather 
ship  data  show  even  less  error,  with  a  slight  summer  bulge  in  the  plot 
likely  due  to  a  warmer  marine  layer  with  a  higher  absolute  water  load. 

We  see  in  these  statistics  the  apparent  effect  of  refraction  gradient  as 
revealed  by  angle  error  calculations  for  different  geographic  locations. 

The  conditions  diffei'  in  the  important  parameter  of  water  vapor  gradient 
(strength  of  the  inversion).  Thus,  there  is  revealed  here  an  additional 
parameter  for  any  predictor  investigation. 


LARGE-SCALE  METEOROLOGICAL.  STUDIES 


In  his  study  of  heat  and  moisture  balance  in  the  marine  layer  Riehl  (1951) 
showed  that  the  net  heat  export  via  the  trade  wind  flow  was  in  the  form  of 
latent  heat.  He  also  found  that  radiative  cooling  at  the  top  of  the  marine  layer 
increased  over  the  average  clear  sky  rate  of  cooling  by  about  33%  in  the  pre¬ 
sence  of  a  stratus  deck.  The  normal  heat  flux  of  the  marine  layer  leaves  a 
residual  of  about  1.  im  cal/day  available  for  net  export  (can  be  monitored  by 
temperature  differential  at  inversion  level).  This  is  about  20%  of  the  heat 
transferred  from  ocean  to  atmosphere  and  is  available  for  modifying  the 
characteristics  of  the  subsiding  air.  For  this  reason  the  amount  of  latent  heat 
surplus  (monitored  via  heat  flux)  could  be  a  parameter  of  interest  in  any 
investigation  of  inversion  activity. 

Of  further  thermal  significance  absorption  of  the  incoming  insolation  at  any 
level  is  proportional  to  the  water  vapor  content  at  that  level,  amd  inversely 
proportional  to  the  total  moisture  content  about  that  level.  Thi  s  places 
maximum  absorption  near  the  base-  of  the  inversion  at  about  the  800  mb  height 
and  suggests  that  the  resulting  thermal  expsmsion  may  contribute  somewhat  to 
those  geometric  changes  of  the  inversion  layer  about  which  we  are  concerned. 
We  therefore  add  water  vapor  at  the  base  of  the  inversion  as  a  parameter  to  be 
monitored  (as  well  as  total  water  vapor  profile  in  the  marine  layer). 

The  principal  output  of  Riehl' s  work,  other  than  explaining  the  changing  geo¬ 
metry  of  the  inversion  layer  in  terms  of  diffusion  of  the  marine  layer  into  the 
dry  air  above,  was  his  explanation  of  the  role  of  vertical  momentum  produced 
by  convection.  Riehl  shows  logically  that  low-level  flow  from  the  east- 
northeast  is  produced  largely  by  the  pressure  gradient.  However,  this  flow 


87 


. . . . 


MONTGOMERY  FIELD,  SAN  01  EGO,  00002.  9, 
AFTER  SHEOOY,  l»G6 


JAN  FEt  MAR  APR  MAY  JUN  JUL  AUC  $EP  OCT  NOV  DEC 


Figure  8-7.  HILO,  HAWAII,  OOOOZ.  0.s  Omr. 
AFTER  SHEDOY,  1966 


g: 


T 


is  checked  by  the  generation  of  vertical  momentum  in  the  convective  cells  of 
the  marine  layer.  The  action  of  this  checking  is  thought  to  be  largely  a  xuiieliou 
of  form  drag  of  the  cells.  This  is  suggested  by  the  relationship  of  wind  speed 
at  various  levels  to  ceil  aivse  at  various  levels.  From  the  standpoint  of 
predicting  inversion  layer  activity,  the  vertiesd  momentum  generated  is 
instrumental  in  the  penetration  of  convective  cells  through  the  stabilized 
inversion  and  thereby  contributes  to  the  height  alteration,  and/or  to  the  total 
destruction  of  the  inversion.  Also,  while  monitoring  the  intensity  of  this 
mixing  at  the  inversion  Interface  through  the  momentum  parameter,  we  may 
wish  to  monitor  the  tramsfer  of  moisture  which  Riehl  shows  Is  not  necessarily 
proportional  to  momentum.  We  therefore  will  wish  to  look  closely  at  the 
Austausch  coefficients  of  momentum  an  i  water  vapor  as  possible  predictors 
of  inversion  layer  activity. 

We  have  briefly  examined  the  work  of  Riehl,  et  al,  and  find  additional 
parameters  tliat  could  logically  be  exercising  control  on,  or  partial  control 
on  the  activity  of  the  inversion  layer.  A  tabulation  of  these  parameters 
plus  those  suggested  by  the  work  of  other  investigators,  will  be  found  in  the 
concluding  secions  of  this  writing.  We  next  examine  the  implications  of 
Neiburger's  mapped  parameter  distributions. 

Perturbations  of  the  inversion  layer  must  be  considered  at  three  general 
scales  of  motion;  macroscale,  or  motions  comparable  to  those  of  large 
cyclonic/anti-cyclonic  weather  systems;  mesoscale,  comparable  to  the  flow 
associated  with  land  and  sea  breezes;  and  microscale,  which  in  this  case  we 
will  view  as  motions  in  the  fine  structure  of  the  inversion  layer  such  as  that 
produced  by  single  cell  convective  penetration  or  gravity  wave  undulations 
found  at  the  marine  layer  and  subsiding  air  interface. 

Turning  here  to  the  macroscale,  we  have  noted  that  the  semi -permanent 
anti -cyclonic  high  (Figure  B-7)  which  sits  over  the  Eastern  Pacific  is  a 
principal  instrument  in  the  creation  of  the  inversion  layer.  One  might  view 
this  high  as  a  large  spiral  descending  vortex.  As  the  air  descends  it  diverges 
outward  and  anti-cyclonlcally  around  the  center  of  the  high.  The  other  input 
necessary  for  the  formation  of  the  inversion  layer  is  the  marine  layer,  which 
is  manufactured  by  moisture  from  the  ocean  surface.  Bearing  in  mind  the 
outward  divergent  flow,  we  note  in  Figure  B-8  the  coincidence  of  the  inversion 
layer  presence  with  the  region  of  this  strong  divergence,  which  suggests 
connection  between  the  two.  Therefore,  we  may  wish  to  add  divergence  to  the 
list  of  possible  statistical  prediction  parameters.  Note  that  the  shaded  area 
of  Figure  B-8  is  coincident  with  the  region  of  strongest  divergence  shown  in 
Figure  1.  The  lessening  inversion  layer  presence  to  the  westward  of  the 
Hawaiian  Islands  area  is  due  to  converging  southwesterly  flow.  The  dropoff, 
of  observed  inversion  presence  along  the  coast  of  Baja,  California  is  due  to 
frequent  north  westward  migration  of  tropical  cyclones  (convergence)  and  the 
dropoff  to  the  north,  off  the  coast  of  Washington,  is  due  to  the  frequent  passage 
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of  winter  storms  (convergence).  South  toward  the  Eqiiator,  the  map  shows 
sn  absence  of  datSr  However,  we  know  from  the  late  1940  cruises  the 
Horizon  and  crest  that  the  inversion  layer  slopes  upward  as  one  proceeds 
south,  and  as  one  approaches  the  zone  of  inter-tropical  convergence  the 
inversion  is  ultimately  destroyed.  The  zone  of  inier-tropieal  convergence  is 
where  the  northeasterly  flowing  trades  of  the  northern  hemisphere  cemverge 
with  the  southeasterly  flowing  trades  of  the  southern  hemisphere*  At  this 
point  of  convergence  the  air  rises  vigorously  and  the  inversion  layer  is 
destroyed  in  the  process. 

Let  us  deviate  for  a  moment  to  explain  the  meteorological  concept  of 
convergence  and  divergence.  Assuming  a  volume  to  remain  constant  as 
we  deform  the  geometry  of  a  cube  of  air,  it  is  clear  from  the  conservation  of 
mass  that  if  we  reduce  the  height  of  the  volume  (vertical  subsidence),  we 
must  expand  laterally  in  either  the  x  or  y  direction  or  both.  Conversely,  if 
we  increase  the  height  of  the  volume,  an  inward  flow,  or  convergence,  must 
take  place.  The  area  of  strong  divergence  in  our  Pacific  anti- cyclone  has 
strong  subsidence  and  outward  lateral  flow.  (Subsidence  is  also  directly 
related  to  vorticity,  or  intensity  of  rotation  of  the  total  vortex. ) 

If  this  subsidence  is  opposed  with  upward  flow  produced  by  the  warm  ocean 
surface  heating  the  air  and  generating  upward  momentum  within  this  convective 
layer,  an  interface  is  formed  between  the  upward  rising  marine  air  sind  the 
downward  subsiding  adiabatically  warmed  air  being  yet  warmer  than  the  ocean 
warmed  air  which  has  cooled  slightly  with  altitude.  Not  only  does  the  temper¬ 
ature  profile  reverse  itself,  but  wbarp  changes  in  water  content  also  occur. 

There  is  strong  indication  from  the  distribution  of  inversion  layer  presence, 
amd  distribution  of  divergence  (or  convergence)  that  a  strong  cause  and  effect 
relation  may  exist  between  these  phenomena.  We  submit  divergence  (or 
subsidence,  or  vorticity:  not  siatiatically  independent  as  a  paranjeter  for 
further  investigation.  What  affects  the  inversion  on  an  eastern  Pacific 
macroscale,  can  also  affect  It  on  a  Vandenberg  local  scale.  Let  us  continue 
with  the  Neiburger  mappings. 

Figure  B-9  is  a  map  of  potential  temperature  diiitribution.  Potential 
temperature  is  a  term  employed  by  meteorologists  to  exclude  compressive 
heating  when  comparing  tow  parcels  of  ajr  from  separate  altitudes.  We 
Note  in  Figure  B-9  that  the  gradient  of  potential  temperature  follows  closely 
the  gradient  distribution  of  ocean  surface  temperature  (Figure  B-10). 
Apparently,  the  marine  layer  temperature  is  being  controlled  specifically  by 
the  temperature  of  the  unuerlyiiig  o(^ean  surface,  at  least  on  the  macro¬ 
scale.  The  cold  sink  off  the  California  coast  is  due  to  well-known  upwelling 
conditions  wherein  cold  ocean  bottom  water  is  brought  to  the  surface.  Note 
in  Figure  B-11  how  the  height  of  the  inversion  base  has  a  distribution  very 
similar  to  that  of  potential  temperature.  Similarities  of  these  distributions 
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suggest  that  a  strong  cause  and  effect  relationship  may  exist  between  these 
paraiiixctcro.  V?c  tliercroi'e  isolale  putvnLIttl  temperature  as  anuiher  param¬ 
eter  to  be  investigated.  Note  also  that  the  inversion  base  slopes  upward  to  the 
wost  from  ahovit,  400  nts  near  the  California  coast  to  more  than  2;  000  m.  near 
the  vicinity  of  the  Hawaiiian  Islands.  The  slope  is  the  greatest  over  the  eastern 
portion  of  the  ocean,  namely  about  1:1000.  As  we  observed  from  sdready 
initiated  field  work,  the  slope  becomes  much  steeper  (1:50)  on  the  Vandenberg 
local  scale.  The  significance  of  slope  is  that  it  alters  the  grazing  angle  of  a 
given  radar  transmission  path  (ray)  as  it  goes  through  the  inversion  layer. 

This  is,  the  ray  may  be  reflected  or  refracted  at  shallow  angles,  depending 
on  the  geometry  of  layer  height  and  statiCKi  height. 

We  liave  excluded  a  diagram  of  the  aver  age  height  of  the  inversion  top  during 
summer.  It  conforms  closely  to  the  topography  of  the  inversion  base, 
suggesting  that  the  dominating  influence  in  determining  these  average  heights 
is  the  same  for  each  case.  However,  an  observed  point  which  does  not  show  in 
the  averages  is  that  a  rise  in  inversion  base  is  often  proceeded  by  a  rise  in 
inversion  top.  We  therefore  add  rate  of  change  of  Inversion  top  to  our  list 
as  a  possible  predictor  of  change  in  inversion  base. 

Figure  B-12  shows  the  distribution  of  inversion  layer  thickness,  with 
minimum  thiclmesa  directly  off  the  coast  of  California  in  the  Vandenberg 
area.  The  area  of  smallest  thickness  and  greatest  temperature  increase  going 
through  the  inversion  layer  reflects  strong  subsidence.  The  greatest 
temperature  at  the  inversion  top  is  in  the  southeast  where  the  inversion  is 
lowest.  A  minimum  sea  surface  temperature  in  the  coastal  area  produces 
cooling  and  settling  of  the  surface  layers,  which  enhances  the  subsidence  and 
divergence  produced  by  the  anti- cyclonic  high  aloft.  This  is  in  contrast  to 
other  regions  where  warming  of  the  sea  surface  produces  strong  cells  of 
warmed  rising  air  which  diminish  the  inversion  contrast.  Where  we  have 
the  strongest  gradients  we  may  also  expect  the  strongest  perturbations. 

This  re-emphasizes  our  need  to  monitor  temperature  gradient  through  the 
inversion  layer  (already  suggested  as  an  indicator  of  heat  flux). 

Figure  E-13  presents  a  summary  of  the  macroscale  cross-section  of  the 
inversion  layer  between  the  Pacific  Coast  and  Honolulu.  First,  note  the 
rapid  rise  of  the  inversion  height  near  the  coast  (excludes  mesoscale 
coastal  upwarp),  with  a  gradual  lessening  of  slope  to  the  westward.  Note 
also  the  liigh  temperature  at  the  top  of  the  inversion  near  the  coast,  pre¬ 
sumably  from  the  strong  subsidence;  and  thirdly,  note  the  gradual  warming 
of  the  entire  air  column  below  the  inversion  as  one  proceeds  westward 
over  an  increasingly  warmer  water  surface. 

I  i 

i  i' 

I  t  We  noted  in  looking  at  this  early  Neiburger  work  that  there  seemed  no  direct 

I  j  relationship  between  windspeed  and  inversion  1  eight.  Other  sources, 

[  I  however.  (Kauper,  1961)  point  out  that  gradien:  winds  in  the  local  Los  Angeles 


95 


. . .  'rn  .  ~  ii  mil  i  i 

. .  '  . . .  ''^Jiii|||iiiii||!|i||||pit|ff  "rip_iii(|inii|ii(||ii|i'!'piii'i[iiii|iit'niiiiiiiiiriiiiliiii|j^ 


IRE  B-12.  AVERAGE  THICKNESS  OF  INVERSION  DURiNG  SUMMER  { HUNDRED'S  OF 
METERS).  AFTER  NEI8UR6ER,  I96S. 


98 


Not  statistically  independent  of  pc*5ntial  ter.ip,  but  could  be  a  better  representation. 


basin  often  orovide  clues  to  the  ehlftlng  height  ef  the  inversion.  Thsrsfors 
we  wish  to  retain  surface  wind  as  a  parameter  for  monitoring. 

SUMMARY  AND  SUGGESTIONS 

We  have  extracted  from  the  literature  the  pertinent  mechanisms  likely  to 
be  controlling  inversion  layer  activity.  The  chief  variables  responsible  for 
the  inversion  layer's  characteristics  are  subsidence  and  resultant  adiabatic 
heating  around  the  anti- cyclonic  vortex  of  the  northeast  Pacific,  heat,  moist¬ 
ure,  and  vertical  momentum  in  the  marine  layer  generated  by  the  sea  surface 
temperature.  Direct  or  indirect  manifestations  of  these  variables  appear 
in  Table  1,  wherein  all  variables  thought  to  have  good  physical  reason  for 
influencing  inversion  layer  activity  are  listed. 

In  suggesting  the  statistical  technique  to  be  employed  in  the  evaluation  of 
these  parameters  as  predictors,  a  form  of  multiple  discriminate  analysis 
should  be  considered  (Enger,  1965).  Computer  analysis  of  single  parameter 
correlation  with  inversion  layer  activity,  and  correlations  of  parameters 
known  to  act  in  concert  should  be  attempted.  The  concept  behind  such 
analyses  is  to  determine  the  discriminatory  power  of  predictors,  evaluated 
by  spatial  separation  as  might  be  envisioned  on  a  scatter  diagram. 

The  variable  to  be  predicted  (e,  g.  height  of  inversion)  is  classified  into 
several  mutually  exclusive  groups.  Drawing  on  the  total  collection  of  inversion 
height- change  data  we  can  use  for  respective  groups  such  distinctions  ais  no 
change  =  0  +  5%;  average  change  =  5%  to  1  standard  deviation;  significant 
change  =  1  standard  deviation.  Expressed  in  symbols  we  might  use  5,  S 
and  >,  respectively.  In  its  most  simple  form  we  can  plot  the  observed  value 
of  our  predictor,  X,  —  say,  change  in  potential  temperature  which  pro¬ 
ceeded  the  inversion  activity,  for  each  case  of  inversion  category  mentioned 
above,  5,  S,  > .  The  average  predictor  value  for  respective  categori^,  Xg, 
X„,  X  ,  might  appear  as  below,  as  would  the  average  for  all  cases,  X. 

5  5  5  555  S  55  SSS  SS  S  S  S  >  B>  > »  >  »  »  >  > 

X.  X.  X  X^ 

One  measure  of  toe  discriminatory  power  of  X  is  the  distcince  between  group 
means,  Xg,  Xg,  X  .  The  greater  the  distance,  the  greater  the  lower. 
Moreover,  the  sp^ad  of  data  points  about  their  respective  group  means  is  a 
second  measure  of  the  effectiveness  of  X  as  a  discriminator.  For  a  given 
distance  between  menas,  the  smaller  the  spread  the  less  the  overlap  would 
be  between  the  groups  and  the  better  the  discrimination.  Spread  is  measured 
by  the  standard  deviations  of  the  respective  groups,  c^,  Og*  o  ^  .  Dis¬ 
crimination  is  maximized  by  the  distance  between  groups  and  minimized  by 
the  spread  of  the  points  about  the  means,  expressed  as 
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\  = 

or, 


rne^Btire  of  distance  aroong  group  meams 
measure  of  spread  about  means 


where  g  represents  5,  S,  >  and  n 


g 


is  the  number  of  cases  in  each  group . 


There  are  refined  techniques  for  selecting  a  minimum  set  of  efficient  predic¬ 
tors  from  a  larger  set  of  possible  predictors,  referred  to  as  screening 
discriminant  analysis.  Many  variations  and  refinements  to  this  general  theme 
are  available,  and  could  be  decided  upon  if  and  when  the  next  phase  of  work  wer 
initiated.  Once  identified,  time  plots  of  strong  parameters  could  be  analy  zed 
for  respective  partial  contributions  to  the  fluctuations  found,  and  time  lag 
between  apparent  cause  and  effect  could  be  similarly  investigated. 


Assuming  controlling  parameters  could  be  statistically  corfiimed,  the  stage 
should  be  well  set  for  mathematical  modeling.  An  excellent  point  of  depai  - 
ture  exists  for  this  phase  of  work.  We  cite  the  healed  island  r  lOdel  of  Lav  ii( 
(1968),  wherein  predictions  of  similar  inversion  activity  are  n.ade  for  pre¬ 
dicting  lines  of  cumulus  build-up. 


There  is  every  reason  to  believe  that  significant  strides  could  be  made  in 
our  ability  to  understand  the  activity  of  the  Pacific  inversion.  A  program  of 
this  nature  would  be  free  of  field  support  requirements.  Access  to  rouUre 
measurements  would  be  the  source  of  necessary  inputs. 


We  point  out  also  that  the  present  work  has  revealed  the  considerable  cor¬ 
relation  between  theoretical  flow  dynamics  and  the  behavior  of  inversion 
distortions  at  Western  Test  Range.  A  Froude  number  approach  would  be 
excellent  for  prediction  and  monitoring  of  local  inversion  anomalies. 
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Anomalous  atmospheric  conditions  at  Western  Test  Range  were  measured  by  aircraft. 
The  inversion  layer  was  found  to  possess  upward  warp  as  it  approached  the  coast,  as 
well  as  dilation,  lensing,  sandwiching,  and  waving.  Flow  geometry,  when  associated 
with  topographic  obstacles,  was  found  to  be  markedly  similar  to  that  produced  on  lab¬ 
oratory  models.  Froude  number  concepts  classified  the  anomalies  accurately.  Re¬ 
fraction  corrections  produced  by  three  possible  data  sources  were  compared  to  assesi 
the  effects  of  the  above  described  anomalies.  When  data  from  full  cross  sections 
along  the  sighting  path  were  employed,  along  with  steering  equations,  elevation  and 
bending  corrections  were  considerably  greater  (2x)  than  those  produced  by  single  poin 
radiosondes  or  aircraft  spiral  measurements.  Moreover,  corrections  produced  by 
spiral  data  increased  to  seaward  where  the  inversion  became  sharply  defined.  The 
validity  of  full  cross  sections  defining  the  propagation  path  was  best  established  by 
examining  ray  paths  for  antenna  elevation  angles  existing  at  the  time  of  signal  loss 
when  a  target  simulation  aircraft  encountered  a  radio  hole.  Ideally,  a  trapped  ray 
could  oe  expected  to  bend  back  to  earth  at  the  aircrafts'  position  in  time  and  space 
when  any  one  of  three  tracking  radars  lost  lock  on  the  target.  A  cross -sectio.a  corr¬ 
ected  propagation  path  defined  this  target  position  better  than  any  other  means  of  de¬ 
fining  the  path.  Strength  and  distinctive  geometry  of  the  inversion  at  the  position  of 
loss  are  displayed.  Mathematics  employed  for  real-time  corrections  at  any  angle 
(including  negative)  are  presented,  as  are  mathematics  for  steering  corrections, 
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